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FOREWORD 


This  final  report  summarizes  the  work  done  by  Raytheon 
Missile  Systems  Division  for  the  TAGSEA  Program  under  General 
Dyneimics  PO  #304490-PB,  prime  contract  no.  N00017-73-C-2244 . 

It  is  submitted  in  compliance  with  Data  Item  A015  and  is 
organized  into  four  volumes  to  ease  handling  and  for  the  conven- 
ience of  the  readers. 

Volume  I,  Clutter  Models,  reports  the  essence  of  the  work 
and  contains  the  models  themselves  which  were  the  prime  objective 
of  the  clutter  portion  of  the  TAGSEA  program;  it  can  be  read  on  a 
stand-alone  basis.  Enough  peripheral  material  is  also  included 
to  provide  a framework  for  a good  understanding  of  the  models. 
Volume  II,  Procedures  and  Output  Forms,  provides  details  and 
explanations  on  methodology  including  the  form  of  the  outputs 
and  the  structures  of  the  clutter  simulation  effort.  Volume  III, 
Supportive  Analyses  and  Outputs,  provides  analytical  back-up  and 
a more  complete  detailed  view  of  the  simulation  software. 

Volume  IV,  Standard  Clutter  Analysis  Outputs,  is  a compilation  in 
various  forms  of  the  mass  of  data  analyzed  during  the  program. 

Each  volume  has  its  own  table  of  contents  which  serves  to  outline 
the  specific  material  presented  therein. 

Raytheon  wishes  to  acknowledge  the  valuable  aid  and  support 
given  by  members  of  the  team  composed  of  personnel  from  NAVSEA, 
APL/JHU,  Technology  Service  Corporation  and  General  Dyn’amics. 

Many  helpful  suggestions  were  made  during  a series  of  critiques 
and  reviews  which  most  assuredly  contributed  to  a better  resul- 
tant output.  The  assistance  received  ranged  all  the  way  from 
general  support  and  overall  guidance  to  specific  supportive 
analyses,  detailed  unpublished  comparative  data,  and  suggestions 
of  exact  forms  of  clutter  models  and  plots  which  would  be  most 
informative  to  the  community  at  large. 
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APPENDIX  A 
SUPPORTIVE  ANALYSIS 


During  the  course  of  the  TAGSEA  program  many  analyses 
were  done  which  were  concerned  with  the  technical  conduct  of 
the  work  rather  than  with  analysis  of  the  resultant  data.  Such 
analyses  dealt  with  the  requirements  for  various  facets  of  the 
programs  including  the  design  of  the  experiments,  application 
of  the  radar  range  equation,  coordinate  transformations,  and 
weighting  matrix  determination  which  are  presented  in  this 
appendix. 

1 . 1 Design  of  the  Experiment 

This  section  discusses  the  considerations  involved  in 
designing  the  clutter-gathering  experiment. 

1.1.1  Objective  of  the  Clutter  Study 

It  was  the  primary  objective  of  the  clutter 
study  effort  to  characterize  sea  clutter  with  emphasis  on  deter- 
mining the  distribution  in  space  and  time  of  the  relatively  rare 
large  radar  cross  section  returns.  Such  a characterization  is 
valuable  for  use  in  future  surface  missile  systems  for  speci- 
fication of  performance  requirements,  design  of  the  missile 
seeker,  as  well  as  simulation  and  evaluation  of  candidate 
seekers.  To  achieve  this  objective  it  was  desired  to  gather 
radar  clutter  data  under  a variety  of  conditions,  to  reduce  and 
analyze  the  data,  and  to  develop  analytical  and  simulation 
models  of  clutter. 
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1.1.2  Data  Requirements 

Sea  clutter  has  been  shown  to  be  representative 
of  a non-stationary  process.  That  is,  statistics  of  the  distri- 
bution function,  such  as  the  mean,  vary  with  time  more  than  can 
be  attributed  to  chance  alone.  It  has  been  suggested  that  the 
slowly  varying  nature  of  the  distribution  can  be  attributed  to 
change  in  the  shape  of  the  coarse  structure  of  ocean  surface. 

A typical  seeker  looking  down  at  the  sea  can  resolve  a patch  on 
the  surface  of  the  ocean  approximately  100  feet  square.  It  is 
known  that  small  patches  exhibit  time-varying  distributions 
with  decorrelation  times  varying  over  a couple  of  orders  of 
magnitude,  depending  on  conditions.  On  the  other  hand,  larger 
patches  can  exhibit  exponential  cross  section  distributions  and 
somewhat  shorter  decorrelation  times.  The  nature  of  clutter 
statistics  are  needed  for  patches  of  sea  100  feet  square  viewed 
at  grazing  angles  between  10®  and  50®.  The  questions  to  be 
answered  are: 

1)  How  often  does  a large  cross  sections  occur? 

2)  How  long  does  it  persist? 

3)  What  is  the  nature  of  rapid  fluctuation  when 
the  mean  cross  section  is  high? 

The  answers  to  these  questions  influence  the 
final  choice  of  parameters  in  clutter  blanking  and  CFAR  circuits 
of  future  seekers. 

The  amount  of  clutter  data  required  to  give 
answers  is  determined  by  the  required  false  alarm  probability. 
Some  typical  seeker  parameters  are  used  for  this  estimate  as 
follows.  If  1%  of  a 4 second  search  period  is  to  be  devoted 
(on  the  average)  to  dismissing  false  alarms,  then  the  required 
false  alarm  probability  per  bin  is  approximately  10  ® for  a 
64  X 64  range-doppler  matrix.  Only  two  out  of  64  doppler  bins 
will  be  likely  to  cause  false  alarms,  so  a false  alarm  prob- 
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ability  of  3 X lO”^  due  to  clutter  would  be  reasonable.  The 
distribution  curve  should  be  measured  up  to  levels  corresponding 
to  this  false  alarm  probability  and  would  require  that  the  level 
be  exceeded  about  16  times  in  order  to  estimate  the  level  to  a 
la  value  of  O.ldB.  This  would  require  16/3  X lO'*’^  = 500,000 
decorrelation  bins  to  be  observed.  The  decorrelation  time  is 
expected  to  be  less  than  0.2  seconds,  so  100,000  bin-seconds 
of  information  would  be  adequate  for  each  condition. 

1.1.3  Approach 

Several  data  gathering  flights  were  made,  in 
which  runs  were  performed  under  the  9 combinations  of  the 
following  conditions: 

Wind  Direction:  a)  Upwind 

b)  Downwind 

c)  Crosswind 

Grazing  Angle:  a)  Low  (11*-15®) .Altitude  1100  ft 

b)  Medium  (22®-32“) , Altitude  2200  ft 

c)  High  (34®-52*),  Altitude  3300  ft 

Each  run  lasted  approximately  five  minutes  and  produced  data  in 
16  contiguous  100  foot  range  gates  using  a side-looking  horn 
as  an  antenna.  (Contiguous  gates  were  chosen  to  provide  a 
map  of  a continuous  patch  of  sea  on  each  run,  and  to  minimize 
the  required  receiver  dynamic  range) . The  output  corresponding 
to  each  range  gate  was  box-carred,  filtered,  and  recorded  on 
tape  using  frequency  multiplexing. 

1.1.4  Data  Gathering 

It  was  desirable  to  gather  data  quickly  in  order 
to  reduce  flight  time.  This  was  achieved  by  viewing  as  large  a 
surface  area  as  possible  consistent  with  aircraft  velocity  and 
radar  average  power.  It  was  also  desirable  to  use  a fixed 
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antenna  thereby  putting  minimum  demands  on  attitude-stability 
of  the  aircraft  since  the  high  sea  state  operations  took  place 
in  rough  weather.  These  considerations  dictated  the  use  of  a 
horn  with  about  60“  beamwidth  and  lOdB  gain.  The  vertical  beam- 
width  covered  the  40®  range  of  grazing  angles  allowing  the  horn 
to  be  mounted  with  a depression  angle  of  30®.  The  horizontal 
beamwidth  permitted  the  aircraft  to  fly  at  a speed  of  as  much  as 
300  knots  (500  ft/sec)  and  still  keep  aliased  spectra  about  lOdB 
below  the  desired  spectrum  in  the  region  of  interest. 


follows ; 


The  data  gathering  system  parameters  were  as 


a)  Radar 


Frequency 

X-Band 

Antenna  Gain 

lldB 

Peak  Power 

lOOOW  at  the 

1 antenna 

Pulse  Width 

0 . 2 Visec 

PRF 

19,320pps 

Noise  Figure 

lldB  at  the 

antenna 

No.  of  Range 
gates 

16 

Range  Gate 

16  parallel 

analog  channels 

Output 

box-carred, 

filtered,  and 

offset  in  doppler  one-quarter 
of  the  PRF 


Range  Gate 
settings 

b)  Recorder 

No.  of  Channels 

Speed 

Freq.  Response 

Record  Time/Reel 

Clutter  Data 

recording 

method 


4500  to  6100  feet 
7 

60  ips 

400Hz  - 750KHZ 
30  min. 

4 frequency-multiplexed  FM 
channels  on  each  of  4 tracks 
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c)  Aircraft 
Type 


1.1.5 


Velocity 

required 

Performance 


250  knots  (420  ft/sec) 


As  the  detailed  calculations  presented  in  Sub- 
section 1.2.1  will  show,  the  received  power  in  a range-doppler 
cell  computes  to  be  -110.4dBm  and  the  noise  power  in  a 5MHz 
bandwidth  is  -96dBm  for  an  lldB  noise  figure.  Thus,  the  clutter- 
to-noise  ratio  is  -110.4- (-96)  = -14.4dB  per  pulse  per  element. 
The  integration  gain  is  128,  or  21dB.  Therefore  the  integrated 
C/N  per  bin  is  21-14.4  = 6.6dB  at  4600  feet  for  = -27.5dB. 

This  is  adequate  for  establishing  the  mean  value  of  a with  less 
than  IdB  error  due  to  noise.  A more  exact  determination  of  the 
statistical  limits  noise  estimation  appears  in  Section  1.3  of 
this  appendix. 

1.2  The  Radar  Range  Equation  and  Resolution 


1.2.1  The  Radar  Range  Equation  for  TAGSEA  - 
Power  Received  in  a Range-Doppler  Cell 


P = 
r 


The  basic  equation  for  power  received  is: 


^ ^t 
4tt  R^ 


1.2-1 


4ir  R^ 


= Transmitted  Power 

= Transmit  Antenna  Gain 

= Range  to  Target  or  Clutter  Point 
= Radar  Cross  Section 
= Receive  Antenna  Aperture  Area 
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Because  a single  X-Band  horn  is  used,  the  gain  on 
transmit  and  receive  is  identical  and  is  related  to  aperture 
area  by  = 4irA^/A^  = G^;  where  A = 0.1  ft  wavelength  at  X-Band. 

Radar  cross  section  o is  the  product  of  reflec- 
tivity a°  and  the  area  A A of  a range-doppler  cell;  the  size  of 
which  is  determined  by  signal  processing.  Microwave  losses  L 
are  about  2dB  each  way  in  the  TAGSEA  radar  hardware.  With  these 
modifications  the  equation  for  power  reduces  to: 


P 


r 


(G/L)^  A^  0°  (A  A) 
(4tt)^  r'* 


1.2-2 


The  Range-Doppler  resolution  cell  size  A A depends 
on  geomstry  as  well  as  signal  processing  and  is  closely  approxi- 
mated hy  the  equation: 

(A/2V)  R (AR)  {Af  ) 

AA  = — - ■■  ^ 1.2-3 

yi-(H/R)^  - (Af^/2V)^ 

' I 

derived  in  Subsection  1.2.2.  Range  resolution  AR  and  doppler 
resolution  A f^  are  fixed  by  signal  processing.  Altitude  H, 

Range  R and  velocity  V are  a function  of  flight  test  conditions. 
Cell  size  increases  slightly  with  altitude  and  doppler  frequency; 
as  indicated  by  the  denominator  terms  in  equation  1.2-3. 

The  antenna  patter  G(a,6)  had  to  be  related  to 
flight  test  geometry  in  order  to  compute  spectral  shape.  The 
expression  for  antenna  gain  (derived  in  Subsection  1.2.3)  was 
modeled  as  being  parabolic  in  decibels. 

Where  G = 9.8dB,  k = 0.00461  deg"^  and  k . = 0.0047  deg"^. 
o az  ’el  ’ 

Angles  (a,g)  are  in  measures  of  degrees.  Combining  equations 
1.2-2,  1.2-3  and  1.2-4  yields  the  result: 
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’zsBcr. 


G^(g,3)  g°(AR)  ^^2-5 

^ {4tt)^  R^(2V)  >/l  “ (H/R)^  -(  Xf^/2V)^ 

Received  clutter  power  in  a range-doppler  cell 
was  computed  in  Table  A-1  for  a nominal  set  of  geometric 
conditions:  R = 4400  ft,  H « 2200  ft,  V = 420  fps  and  zero  doppler. 
For  this  particular  choice  of  parameters  the  depression  angle  is 
30°  and  antenna  gain  is  maximum;  i.e.,  a = 0 = g. 

In  Subsection  1.2.3  it  is  shown  that  the  antenna 
steering  coordinates  (a,g)  are  functionally  related  to  geometry 
by  the  vector  relation; 


tan  a 

0 

sin (6g-0^)cosy 

0 

= cosij; 

1 

= 

cosC 

cos (0g-0^)cosy 

+ sinC 

siny 

where ; 

tan  g 

-sinC 

cos (0g-0^)cosy 

+ cosC 

siny 

1.2-6 


e 

g 

«d 

Y 


= azimuth  angle  of  the  range  vector  in  the  ground 
plane 

= aircraft  drift  angle 

= depression  angle  of  the  range  vector 

= depression  angle  (30°)  of  the  antenna (X-Band  horn) 


tanij^ 


2 

tan  a 


tan^g 


The  coordinate  rotation  angles  are  measured  with  respect 

to  a system  determined  by  the  velocity  vector,  altitude  vector 

and  zero  doppler  frequency.  The  depression  angle  is  related  to 

2 

range  and  altitude  by  siny  = H/R  or  cosy  = 1-(H/R)  . Azimuth 

in  the  ground  plane  is  defined  by  sinO^  = ( Xf^/2V)  cosy  (see 
Subsection  1.2.2). 
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TABLE  A-1 

COMPUTATION  OF  POWER  RECEIVED  IN  A RANGE-DOPPLER  CELL 

P.  G^(a,B)  A^o»(AR)  .) 
p _ t a 

(4ir)^R^(2V)  yi-(H/R)^  - ( Xfy2V)  ^ 

H = 2200  ft,  R = 4400  ft,  a = 3 = 0®,  = 0 

Parameter  + - 

P^  = 1000  watts  60dBm  60dBm 


G^(a,6)=  2 X lldB  @ a=0=6  22 

30  log(0.1  ft)  = -30dB  ft^  30.0 

= -27.5dB  @ 30®  depression,  SS-3  27.5 

AR  = 110  ft  +20.4dB  ft  20.4 

Af^  = 10  log(14KHz/128)  = 20.4dB  20.4 

(4tt)^  30  log(4TT)  = 33dB  33 

R^  30  log(4000)  = 109. 9dB  ft^  109.9 

(2V)  10  log(2X420)  = 29.4dB  ft/pec  29.4 

(H/R)^',  5 log(l-0.5^)  = 0.624dB  0.6 

L = 2dB  4.0 

123.4  233.8 

P^  =123.4  - 233.8  = -110.4dBm 
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Equations  1.2-4,  1.2-5  and  1.2-6  were  combined 
to  plot  the  curves  in  Figures  A-1  (SS-3)  and  A-2  (SS-5)  for 
bounding  conditions  of  range  and  altitude.  These  curves  were 
normalized  to  zero  decibels  at  zero  doppler  and  replotted  in 
Figure  A-3.  Because  drift  angle  was  set  to  zero  the  curves  are 
symmetrical  about  zero  doppler. 

Values  of  used  for  plotting  Figures  A-1  and 
A-2  were  taken  from  data  provided  by  APL  in  Figure  A-4 . 

1.2.2  Range-Doppler  Surface  Geometry 

Consider  the  geometry  in  Figure  A-5. 


Figure  A-5  Flight  Test  Geometry 
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DOPPLER  FREQUENCY  IN  HERTZ 
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d 


DOPPLER  FILTER  NUMBER 


Figure  A-3  Normalized  Power  Received  as  a Function  of  Doppler  Frequency 
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Figure  A-4  Normalized  Backscatter  Coefficient 
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The  range  rate  is  given  by: 

V.R  Vx 

R - V cos  0 = |_|  = R 


1.2-6 


Doppler  frequency  is  related  to  range  rate  by  the  equation. 


- 2R  2V  « 

'd  “ — = T 


1.2-7 


The  above  two  equations  can  be  combined  to  provide  an  expression 
for  isodop  contours 


X ^ _2d 

R 2V 


2 ^ 2 ^ 2 
/X  + y + z 


1.2-8 


The  isodops  intersect  the  clutter  surface  at  altitude  z = H.  If 
a polar  coordinate  system  on  the  clutter  surface  is  chosen  as 
indicated  in  Figure  A-6  the  isodop  equation  becomes 


R sin0 


2 2 
R + H 

g 


= - 


2 

(H/R)‘‘  sine 


1.2-9 


The  isodop  equation  was  used  to  plot  the  ground  maps  in 
Figure  A-6  as  well  as  those  in  the  body  of  the  report. 

The  area  of  the  clutter  patch  in  Figure  A-7  is 


dA  = (R  de  )dR 

g g g 

It  is  easily  shown  that 


1.2-10 


dRg  = d ( V ) = 


1.2-11 


/r2  - 
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Figure  A-6  Surface  Map  Overlays 
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and  from  equation  1.2-9  that; 


de  = d 

g 


. -1  / ^d/2V 

sin  ' 


(H/R) ^ 


1.2-12 


( V2V)  df 


^ 2 2~* 
.1  - (H/R)^  - (Af^/2V)^ 


UNCLASSIFIED 


UNCLASSIFIED 


Substituting  equations  1.2-11  and  1.2-12  into  equation  1.2-10 
provides  the  resultant  expression  for  small  areas  as  a function 
of  range  and  doppler.  Replacing  differentials  by  deltas  yields: 


a A 


2V 


XR  (AR)  (LfJ 
d 

V 1 - (H/R)2  - 


( Af^/2V) 


1.2-13 


Coordinate  rotations  can  be  derived  directly 
from  Figure  A-5.  If  a unit  vector  (010)  along  the  range 
to  a clutter  point  is  to  be  rotated  to  the  coordinate  system  of 
the  velocity  vector,  the  transformation  is: 


COSOg 

-sInOg 

0 


sin0 

g 

COS0 

g 

0 


o 

o 



0 

0 cosy  -siny 

1 

0 siny  cosy 

0 

cosy 

”"®g 

(Xf^/2V) 

u = 

cosy 

COSOg^ 

= 

V 1-(H/R)2-  (Xf^/2V)"  ' 

siny 

H/R 

where:  siny  = H/R 


cosy 


-^1  - (H/R)2' 


sln6 

g 


(Afj/2V) 

I ^ 

V~  (l!/!<) 


COS0 

g 


(Xf^/2V)2 
1 - (H/R)- 
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1.2.3  Antenna  Coordinate  Transformation 

Equations  for  rotation  of  points  in  antenna 
space  to  a horizontal  cartesian  (x,y,z)  coordinate  systems 
based  on  x along  the  velocity,  y along  zero  doppler  and  z along 
the  vertical  are  derived  from  Figure  A-8.  The  unit  vector  u^ 
in  antenna  space  is  transformed  to  velocity  space  by  the 
relation: 


V 

or,  using  the  matrix  notation,  I'] 


V 


C = Antenna  Depression 
Angle  (X-Band  Horn) 

0^  = Aircraft  Drift  Angle 

Clutter 

Patch 


s j nG 


0 1 


cosOj  0 0 cosf;  -sinC  cosi(/  1 1.2-15 


i jjl)  sin£.  cos^ 


tana 


[tan8 


Figure  A-8  Relationship  Between  Antenna  Space  and  Velocity  Space 
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The  matrix  relation  for  transforming  the  range  vector  into  the 
velocity  vector  derived  in  Subsection  1.2.2  is  given  by; 

= [Qgl  [y]  Uj^  1.2-17 

where 

= [0  1 0] 

The  problem  is  to  find  an  expression  for  the 
antenna  vector  in  terms  of  the  range  vector.  Thus,  using  matrix 
algebra  and  the  previous  two  equations  gives  the  result: 

‘'A  ' [OgJ  [Y]u^  1.2-18 

Because  the  0,  matrix  and  0 matrix  are  of  the 
d g 

same  functional  form  it  is  easily  shown  that: 
i®gi  ■ [®g  ■ ®d) 

Therefore,  equation  1.2-18  reduces  to: 

“ [5]'^  [®g  * [Y]Ug  1.2-20 

Expansion  of  equation  1.2-20  yields  the  final 
expression  in  antenna  space: 


10  0 If  cos(0  -0.).  8ln(e  -0.)  Olfl  0 0 1 [o 

u = 0 cosC  8lnC  -8ln(0  -0^)  coa(0  -0^)  0 0 1 

0 -sine  cose  0 * 0 * 10  0 
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The  matrix  products  gives  the  result: 


tana 

0 

sln(9  -0^)co8v 

g d 

— « 
0 

U^-COSl(/ 

1 

m 

cos  C 

cos(e  -0,)cosY  + sln£ 
g 0 

slnY 

tanB 

-Sint 

cos(0  -0,)cOSY  + cost 
” 8 ^ — 

slnY 

1.2.4  Transformation  to  surface  Coordinates 

The  clutter  data  was  gathered  by  a monostatic  radar 
mounted  on  a moving  airplane.  Consequently  the  coordinate 
system  with  the  data  was  that  of  the  aircraft.  Radar  returns 
were  collected  and  sorted  into  bins  by  their  range  and  doppler 
(range  rate) from  the  moving  aircraft. 

Let  us  assume  constant  velocity  and  altitude 
for  the  aircraft.  Now  defining  the  following  aircraft  base 
and  coordinates 

Xa  along  the  velocity  vector  (horizontal) 

Ya  horizontal  and  perpendicular  towards  the  clutter 

collection  side 

Za  altitude, 

the  range  rate  R to  any  stationary  point  from  the  moving  plat- 

o 

form  is  clearly 


R = V cos  a 


1.2-22 


where  cos  a is  the  half  cone  angle  between  range  and  aircraft 
velocity  vector,  see  Figure  A-9. 


R 


2 


The  following  relationships  are  easily  shown. 
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V X 


a 
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R 


R 


1.2-24 
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Aircraft 

Position 


Figure  A-9  Aircraft  Relationship  to  a Surface  Clutter  Cell 


For  any  particular  cell  we  have  coordinates 
R,  R and  are  presuming  known.  Solving  for  the  rectangular 
coordinates  of  the  clutter  with  respect  to  the  aircraft  we  have: 


1.2- 25 

1.2- 26 


If  we  now  assume  constant  velocity  V and  a 
surface  coordinate  system  centered  under  the  aircraft  at  t=0; 
X,  Y,  Z at  time  t are  given  by: 


X = X + Vt  1.2-27 

a 

Y = Y 1.2-28 

a 


Z = Z = altitude 

cl 
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1 
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Our  mapping  of  data  into  surface  coordinates 
is  given  by  equations  1.2-25  through  1.2-29  where  we  measure: 

R = from  range  gate  identification 

R = from  doppler  cell  identification 

t = from  frame  index  ^ 

We  assume  from  pilot  control: 

Z = aircraft  altitude 

V = aircraft  velocity  (surface  velocity) 

1.2.5  Transformation  in  Terms  of  Transferred  Data  Form 
Inputs  available  from  flight  plan  and  data  logs 
are;  V,  aircraft  velocity  and  Z,  aircraft  altitude.  From  the 
data  tape  the  following  inputs  are  available;  A = power  of  return, 
I = range  gate  number,  J = doppler  cell  number,  K = frame  number 
(number  of  FFT's  since  start  of  run). 

The  transformation  is  computed  as  follows: 

A (I,  J,  K) rA  (X,  y,  t) 

t = t^  + K(At) 

R = R^  + 1(AR) 

• * 

R = R^  + J(AR) 

X'  = R R/V 
X = X'  + VT 

Y = Vr^-X'^-Z^ 

AT  is  FFT  frame  time  (data  reduction  parameter) 

AR  is  range  gate  dimension  (system  parameter) 

A R is  doppler  cell  size  in  velocity  unit 
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as : 


Note  that  the  doppler  cell  size  can  be  defined 


A R 


c 


Where ; 

AF  = frequency  cell  size 
f^  = range  frequency 

c = speed  of  light 


1.2.6  Radar  Resolution  Characteristics 

Radar  range  and  doppler  resolution  characteristics 
are  described  in  this  section.  Range  resolution  is  primarily 
determined  by  transmitted  pulsewidth  and  receiver  video  filters 
ahead  of  the  sample  and  hold  circuits.  Doppler  resolution  is  a 
function  of  the  data  reduction  process  and  its  associated  Fast 
Fourier  Transform  (FFT)  signal  processing  equipment. 


The  transmitted  wavefojnn  envelope  is  closely 
approximated  by  a rectangular  CW  pulse  with  width  T = 0.2usec 
time  duration.  The  two-sided  IF  bandwidth  of  the  receiver  is 
6.6  MHz  with  spectral  characteristics  determined  by  a two-pole 
Bessel  filter.  Time-bandwidth  product  for  the  filtered  waveform 
is  therefore  BT  = 6.6MHz  X 0.2usec  = 1.32.  The  effective  pulse- 
width  has  been  computed  as  1.2T  = 0.24usec*. 


* L.  W.  Brooks,  "Equivalent  Pulsewidth  of  the  TAGSEA  Data 
Gathering  System",  Technology  Service  Corporation  memo 
TSC-W2-65  B50711,  dated  27  April  1976. 
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An  ambiguity  diagram  for  the  transmitted  pulse 
doppler  waveform  is  sketched  in  Figure  A-10.  The  PRF  used  for 
transmission  of  the  pulse  doppler  waveform  was  19.32KHz.  Range 
ambiguities  in  the  "bed  of  nails"  are  separated  by  25,472  feet, 
or  4.19  nmi.  Velocity  ambiguities  are  separated  by  974  feet/sec. 
The  darkened  area  in  Figure  A-10  shows  the  range-velocity 
region  mapped  by  the  system.  It  is  extended  over  1600  feet 
in  range  and  176.4  ft/sec  in  velocity.  The  first  range  ambiguity 
is  at  25.5KFT  and  the  first  velocity  ambiguity  is  at  974  ft/sec. 
Therefore,  the  mapped  region  is  clear  from  all  ambiguities  in 
range  and  velocity. 

Figure  A-11  shows  an  enlarged  view  of  the 
ambiguity  diagram  in  which  both  the  waveform  characteristics 
(after  signal  processing)  and  the  mapped  region  are  displayed  in 
greater  detail.  Range  resolution  for  an  effective  pulsewidth 
of  0.24usec  is  118  feet.  Velocity  resolution  calculations 
require  an  understanding  of  the  data  reduction  process. 

Briefly,  pulse  doppler  echoes  were  applied  to 
a sample  and  hold  and  boxcar  integrator  at  a rate  equal  to  the 
radar  PRF  (19.32  KHz).  The  integrated  data  was  recorded  on  an 
analog  tape  for  each  range  gate.  Later,  during  data  reduction 
playback,  each  range  gate  was  played  back  and  re-sampled  at  a 
reduced  rate  of  14  KHz.  Continuous  trains  of  128  samples  were 
processed  without  weighting  by  a Fast  Fourier  Transform  (FFT) 
filter  bank.  Thirty-two  out  of  the  128  filters  output  from  the 
FFT  were  used  for  spectrum  analysis.  The  spacing  between  doppler 
filters  was  14KHz/128  = 109.375  Hz.  Because  the  weighting  was 
uniform,  each  FFT  filter  response  was  essentially  a sin  x/x 
function  with  doppler  resolution  equal  to  0.886  of  the  filter 
spacing,  or  96.57  Hz.  Velocity  resolution  at  X-band  is  given 
by 

''rzS  = I ^d(RES)  = 
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-10  Radar  Ambiguity  Diagram 


Figure  A-11  Enlarged  View  of  the  Ambiguity  Diagram 
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The  range  velocity  waveform  ambiguity  diagram 
is  compared  to  the  mapped  region  in  Figure  A-11.  Range  resolution 
(118  ft)  is  slightly  greater  than  the  range  gate  spacing  (100  ft) . 
The  velocity  resolution  (4.88  ft/sec)  is  slightly  less  than  the 
spacing  of  velocity  cells  (5.5125  ft/sec).  The  mapped  region 
spans  16  cells  or  1600  feet  in  range  and  32  cells  or  176.4  ft/ 
sec  in  velocity. 

1 . 3 Statistics/Confidence 


1.3.1  Statistics  of  Receiver  Noise 

This  section  is  concerned  with  estimating  the 
mean  noise  level  of  the  TAGSEA  radar  system  and  obtaining  con- 
fidence intervals  for  each  estimate.  At  the  end  of  each  run  a 
10  to  15  second  interval  of  a calibration  signal  was  recorded 
and  another  interval  of  10  to  15  seconds  was  noted  with  the 
transmitter  off.  This  second  period  was  an  observation  of  the 
noise  level  of  each  range  gate  in  the  system  during  that  run. 

The  noise  record,  the  calibration  signal  record,  and  the  record 
of  the  setting  of  various  meter  readings  were  all  kept  as 
possible  references  for  determining  the  gain  of  the  receiver 
range  gates.  The  received  level  was  adjusted  before  each  run 
so  the  average  clutter  power  level  fell  between  predetermined 
intervals. 

The  effect  was  to  put  the  noise  rather  than  the 

calibration  signal  more  in  the  dynamic  range  of  the  clutter 

return  in  most  range  gates.  This  made  the  average  noise  power 

a more  stable  reference  than  the  calibration  signal  and  a more 

accurate  measurement  of  the  range  gate  gains;  doing  this  allowed 

a better  calculation  of  a . 

o 
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In  processing  of  the  clutter  data  from  analog 
to  digital  to  form  histograms,  the  signal  strength  was  recorded 
only  as  a value  between  certain  limits  without  regard  to  exactly 
where  in  that  interval  it  fell.  This  "grouping"  of  the  data  has 
little  effect  on  the  uncertainty  of  the  estimate  of  the  average 
noise  level  for  large  noise  levels.  (See  Figure  A-12) . For 
low  noise  levels  the  grouping  has  a more  serious  biasing  effect 
in  the  estimate  of  the  average  noise  level.  The  grouping  of 
the  data  was  performed  by  collecting  the  sample  values  into 
bins  (or  class  intervals)  depending  upon  their  magnitude.  There 
were  1024  class  intervals  and  is  a sample  fell  between  levels 
0 and  1 but  was  not  equal  to  1,  it  was  counted  as  a value  in 
bin  1.  In  an  effort  to  maintain  a constant  average  clutter 
return  from  range  gate  to  range  gate,  the  gains  of  the  closer 
range  gate  filter  were  low  enough  in  most  runs  to  place  most  of 
the  noise  signal  in  bin  1.  See  Figure  A-13.  The  question 
naturally  arises  as  to  the  accuracy  of  any  estimate  of  the 
average  noise  level  for  those  range  gates.  Due  to  the  grouping, 
the  average  power  level  can  be  anywhere  between  0 and  1.  Con- 
fidence about  the  estimate  of  the  mean  of  receiver  noise  is 
established  using  the  concept  of  confidence  intervals. 

The  most  notable  results  coming  from  the  use 
of  confidence  interval  estimation  is  the  relatively  small  error 
that  is  suffered  for  a high  degree  of  confidence.  If  we  have  at 
least  one  value  of  noise  power  above  class  interval  1 we  can 
make  an  estimate  of  its  mean  and  feel  confident  that  99  out  of 
100  of  these  estimates  will  be  in  error  by  at  most  1.6dB.  This 
is  a statement  which  can  be  made  knowing  only  that  we  have  at 
least  one  hit  above  bin  one.  Smaller  errors  naturally  result 
from  knowing  exactly  how  many  hits  above  class  one  there  are. 
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Figure  a-13  Histogram  of  Noise  with  Histogram  Mean  About  0.5 
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Figure  A-12  Histogram  of  Noise  with  Mean  About  Bin  15 
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The  succeeding  analysis  is  based  on  the  assumption 
that  the  receiver  power  is  exponentially  distributed  and  hence 
determination  of  only  one  parameter,  the  mean  of  the  noise  power, 
will  completely  describe  the  distribution. 

1.3.2  Interpretation  of  the  Results  of  Receiver 

No3.se  Analysis 

A complete  understanding  of  how  to  interpret 
the  results  is  needed  before  a summary  of  the  analysis  is  given. 
The  object  is  to  make  some  inference  about  where  a parameter 
of  a one  parameter  distribution  lies,  in  this  case  the  mean  of 
the  distribution.  We  want  to  make  such  an  inference  with  a 
certain  degree  of  confidence,  which  will  be  measured  by  a 
"confidence  coefficient".  The  mean  of  the  noise  during  a 
particular  noise  histogram  is  a fixed  quantity,  not  a random 
variable.  This  precludes  inferences  of  the  kind,  "the  prob- 
ability of  the  mean  being  situated  between  two  values  is  P, 
where  P is  a confidence  coefficient".  A statement  of  that  type 
would  make  sense  only  if  the  mean  were  a random  variable.  How- 
ever, it  is  also  realized  that  from  noise  calibration  to  noise 
calibration  the  receiver  gain  may  and  does  vary  due  to  system 
variations  and  gain  settings.  The  hope  is  to  make  an  inference 
about  the  mean  without  having  to  hypothesize  about  its  random 
and  non-random  nature.  The  concept  of  "confidence  intervals" 
is  one  way  to  handle  this  problem. 

Confidence  intervals  allow  us  to  make  the 
following  statement:  For  a pre-assigned  confidence  coefficient 
y (0<Y<1)  , we  will  count  the  number  of  occurances  x^  above 
interval  1 and  using  y and  x determine  lower  and  upper  limits 
(jj_and  ~)  and  claim  that  the  true  mean,  lies  between  and 

17'.  Either  the  mean  is  in  that  interval  or  it  is  not,  but 
in  the  long  run  our  inferences  about  where  the  mean  lies  will 
be  in  error  by  at  most  100  (l-y)  percent  of  the  time. 
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If  we  choose  y close  to  1,  we  will  be  in  error 
a small  percentage  of  the  time.  The  price  paid  for  a small 
error  is  that  the  confidence  interval  may  then  be  large. 

i 

i 

1.3.3  Conclusions  on  Statistics  and  Confidence  | 

of  Receiver  Noise 

The  results  of  the  confidence  interval  estimation 
is  best  described  by  referring  to  Figure  A-14  which  tabulates 
what  the  confidence  interval  (jj,  y)  is  for  specific  hits  above 
interval  1 and  for  various  confidence  levels.  As  expected, 
the  figure  does  show  that  the  lower  the  confidence  level,  the 
narrower  the  confidence  interval  is.  For  the  values  of  y in 
Figure  A-14  the  largest  error  would  be  committed  if  y = 99%  and 
the  number  of  hits  (rt^ ) equal  1.  In  this  special  case  if  the 
true  mean,  y^,  were  within  the  confidence  interval  .070£y^  £ .127, 
an  estimate  of  0.10139  for  the  mean  is  in  error  by  at  most  1.6dB. 

This  is  the  maximum  error  at  the  99%  confidence  level  or  less 
given  one  hit.  The  estimate,  y* , of  the  mean  is  obtained  by 
the  equation 


ln{  A) 

where  A is  the  ratio  of  hits  above  interval  1 to  total  hits.  If 
for  each  run  we  form  the  ratio  A and  use  the  above  equation  to 
get  an  estimate  y*  of  y^  we  note  that  at  least  99  out  of  100  of 
these  estimates  will  be  in  error  by  at  most  1.6dB  assuming  that 
we  have  at  least  one  hit  above  interval  1. 

Figures  A-15,  A-16  and  A-17  give  contours  from 
which  other  confidence  levels  can  be  calculated  for  specific 
confidence  coefficients.  These  figures  are  for  y = 0.99, 

0.90,  0.6826  and  hits  less  than  or  equal  to  100. 
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Figure  A-16  Plot  of  90%  Confidence  Interval 
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Number  of  Hits  Above  Interval 


Figure  A-17  Plot  of  68.26%  Confidence  Interval 
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Using  Figures  A-15,  A-16  and  A-17,  the  following 
procedure  is  suggested  to  obtain  confidence  intervals.  For  a 
certain  number  of  hits,  x^,  in  the  range  given  in  each  figure, 
draw  a horizontal  line,  Y,  through  the  point  (0,  x^)  in  the  figure 
with  the  desired  confidence  coefficient.  Where  Y crosses  the 
two  contours  ( y,)j  ) and  T2(Yfy)»  drop  down  to  the  horizontal 
axis  and  read  and  ~.  For  example,  with  20  out  of  19,200 
values  above  class  interval  1 and  confidence  at  the  90%  level 
go  to  IigureA-15.  Draw  a horizontal  line  Y through  (0,  20)  and 
from  the  intersection  of  Y and  T2(*90,  p)  drop  down  to  the  y axis 
and  read  0.137  for  j^.  The  upper  confidence  limit  y = 0.154  is 
obtained  by  dropping  down  from  the  intersection  of  Y with 
T^(.90,  y)  . 

If  confidence  intervals  for  hits  greater  than 
those  given  in  Figure  A-15,  A-16  and  A-17  are  desired,  the 
approximations  given  in  Subsection  1.3.6  can  be  used.  See 
equations  1.3-30  and  1.3-31. 


For  a large  number  of  hits  above  interval  1 , 
the  error  introduced  in  estimating  the  mean  using  the  histogram 
mean  is  no  longer  large.  As  will  be  seen  a relationship  between 
the  noise  mean  y and  the  histogram  mean  y^  under  certain  assumptions 
is  given  by 


I 


ln( 


+1. 


2^  -1 


1 A 

-<y<  CO 

1.3-2 


y = 


0 


A 

y = 


1 

2 


The  error  in  estimating  y using  the  above  transformation  for 
large  number  of  hits  above  interval  1 can  more  easily  be  obtained 
by  using  the  fact  that  the  standard  deviation  ^ of  the  sample 

X 

mean  of  the  grouped  data  is 
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1.3-3 


where  N is  the  number  of  data  points  grouped  and 
U is  the  actual  mean  of  the  noise  power. 

1.3.4  Correction  on  the  Average  Noise  Estimate 

Grouping 

To  determine  estimates  of  the  probability  density 
function  for  clutter  returns  the  data  was  grouped  to  form  histo- 
grams. The  noise  signal  was  also  grouped  and  the  mean  of  the 
noise  power  was  estimated  using  the  sample  mean  of  the  noise 
histogram.  The  grouping  was  performed  by  collecting  the  data 
points  into  bins  or  class  intervals  according  to  the  following 
scheme . 

If  X is  a sample  value  such  that  k-1  ^ x < k 
for  1 k 1023  then  x is  counted  as  a value  in  bin  k,  but  for 
x>1023  it  is  counted  as  a value  in  bin  1024.  To  form  histograms 
the  number  of  points  in  each  bin  is  multiplied  times  a represen- 
tative value  of  the  data  in  that  bin  called  the  class  mark.  If 
the  midpoint  of  each  bin  is  chosen  as  the  class  mark  then  for 
bin  k the  class  mark  is  k-|-  . A class  mark  of  the  left  endpoint 
of  each  bin  was  used  as  a class  mark  for  certain  data  reduction 
but  the  following  equations  are  developed  only  for  a class  mark 
at  the  midpoint. 

If  we  assume  the  data  x is  sampled  from  a 
stationary  random  process  with  a continuous  distribution  F(x) 
and  probability  density  function  f (x)  then  grouping  the  data 
implies  we  are  no  longer  sampling  the  original  random  variable  x 
but  one,  y,  ^/hich  is  now  discrete  taking  only  values  of 
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function  (pdf)  of  y be  f (y)  and  assume  y is  not  bounded  by 
1023^  then 

f(y)  = Pu  for  y = , k = 1,2, 3, 4..., 


where 


-L 


) dx 


k = 1,2,3. 


1.2-4 


If  we  take  N independent  sampled  values  of  x, 
say  x^,  . . . , x^  and  group  them,  they  give  N independent  sample 
values  of  y.  Forming  the  sample  mean,  y,  of  the  grouped  values 
gives  an  estimate  of  the  mean  fj  of  y and  y is  a biased  estimate 
of  the  mean,  y,  of  x.  If  we  assume  that  the  noise  power,  x, 
is  exponentially  distributed  then  the  relationship  between  y 
and  y is 


00 

^ /2k-l.  ^ 

P - (—y ) Pr 

k=l  ^ * 


or 


!>=  f t i 

k=l  ^ k-1  ^ 


dx 


This  simplifies  to 


1.3-5 


1.3-6 


y = 2 


1 (1+e^  ) 


(l-e‘-P  ) 


and  solving  for  y gives 
_1 

|ln<4i) 

v = ' 


M 0 


1 A 

2 < P < 


y = 2 
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The  standard  deviation  of  the  sample  mean  y is  easily  seen  to  be 


/ 

/ 1 e ^ 

" V N TT 

y f (1  - e 


1.3-9 


Where  N is  the  number  of  samples, 


1.3.5  Confidence  Interval  Analysis  of  Noise  Mean 
To  obtain  confidence  regarding  the  mean  of 
receiver  noise  from  noise  histograms,  the  theory  of  interval 
estimation  was  employed. To  establish  a "confidence  interval" 
for  the  mean,  some  probalistic  model  of  receiver  noise  is  needed. 
The  exponential  model  for  noise  power  is  known  to  be  a good 
approximation  and  in  the  following  analysis  it  is  assumed  that 
the  noise  voltage  is  Rayliegh  distributed.  This  assumption 
implies  that  determination  of  only  one  parameter,  the  mean  of 
the  noise  will  completely  describe  the  distribution.  This  one 
parameter  distribution  makes  determination  of  confidence  interval 
easier  to  handle  and  lends  itself  to  fairly  standard  statistical 

analysis. (3) 

If  X is  the  random  variable  y describing  noise 
power,  then  the  pdf  (probability  density  function)  f(x)  is 


f(x)  = 


T e 


- X. 

y 


X > 0 


X < 0 


1.3-10 


(1)  H.  Cramer  (1946),  "Mathematical  Methods  of  Statistics", 
Princeton  University  Press. 

(2)  M.  G.  Kendell  and  A.  Stuart,  "The  Advanced  Theory  of 

, Statistics",  Vol.  II,  2nd  Edn,  Chapter  20,  pp.  98  - 120, 

(3)  Samuel  S.  Wilks,  (1963),  "Mathematical  Statistics", 

John  Wilkey  & Sons,  Inc. 
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where  vi  = E (x)  . 

Considering  only  hits  above  interval  one  without  recording  where 
above  one  a hit  falls,  we  obtain  as  the  probability  p of  getting 
a hit 

' I 


'dy 


1.3-11 


Evaluating  the  far  left  side  of  equation  1.3-11  gives 
1 


1.3-12a 


Equation  1.3-12a  is  important  since  it  relates  to  the  mean,  y, 
of  the  exponential  distribution  to  the  probability  p through  an 
invertible  function.  The  task  of  establishing  confidence  inter- 
vals for  the  parameter  p with  a pre-assigned  confidence  level, 
y,  can  be  translated  to  one  of  estalishing  confidence  intervals 
for  the  parameter  p with  the  same  confidence  level  y and  then 
translating  back  by  using  equation  1.3-12a  or  equivalently 


P = - 


ln(p) 


o<p< 


1.3-12b 


Obtaining  an  efficient  estimate  of  p requires 
only  counting  the  number  of  hits  above  1 in  n independent 
samples  x^^,  X2 , 

given  a value  of  x,  let  A (x)  be  the  function. 


. . , x^  and  then  dividing  by  n.  Therefore 


A(x)  - 


1 i f 

0 if  X<1 
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If  X is  a random  variable  with  pdf  given  by  equation  1.3-10,  then 
A (x)  is  a random  variable  with  pdf  given  by 


. I I'  “■  r 

f A (A  ) = \ 1-p  if  A 


1.3-14 


o 3 f A / o , 1 


where  p is  given  by  equation  1.3-11.  The  expected  value  and 
variance  of  A are 


li{A}  «=  p iiiul  0 “ P(l-P) 


1.3-15 


^1'  ^2*  • • • ' ^ independent  samples  of  receiver  noise 

with  each  x having  the  same  pdf  as  x (equation  1.3-10).  Forming 
the  sum  S 

S = A(x^)  +A  (xp 

gives  a random  variable  x which  measures  the  number  of  hits  above 
class  interval  1.  Dividing  S by  n gives  the  sample  mean~  of  the 
independent  samples  A (x^^)  , A (X2)  , . . . , A (x^)  , i.e.  , 


n 

^ ^ I A(X  ) 

" k = l 


1.3-16 


For  particular  sample  values  x^^,  X2,  ...,  is  the  relative 

frequency  of  the  number  of  hits  above  class  interval  1 in  n 
hits.  The  cdf  (cumulative  distribution  function)  ofT  is  just 
the  binominal  distribution  given  by 


r(r;p) 


io(0  " 


(i  pV 
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To  apply  the  theory  of  confidence  intervals  to 
the  pareuneter  p,  let  y,  0<y<l,  be  the  confidence  coefficient  and 
let  0^  = . Let  T^(Y/P)  and  T2(y»p)  be  two  integer  valued 

functions  of  y ^nd  p such  that  for  a fixed  y ^nd  0<p<l  we  have 
the  following  inequalities  using  equation  1.3-17. 


F(Tj^/n;p)  < a and  F(Tj^/n;p)  ^ a 


1.3-18 


for  (Y»P)  and 


l-F(T2/n;p)  < o and  l-F(T2/n;p)  ^ a 1.3-19 

for  T2 (y»P) • 

It  can  be  shown  that  for  the  binominal  distri- 
bution with  parameter  p functions  T, , T«  given  by  equations 

^ ^ 13)  - 

1.3-18  and  1.3-19  do  exist  for  every  value  of  p.  Since  A 

is  a discrete  random  variable,  Tj^(y,p),  T2(y»P)  are  step  function 

with  properties  that 

a)  T2(y»P)>Tj^(Y/P)  for  0<p<l,  0<y<1 

b)  T2»  are  non  decreasing  functions  of  p as 
illustrated  in  Figure  A-18. 


Note  that  for  a fixed  value  of  p,  say  p^,  we 


have 


’r  {■'l  (S)<P, 


o'’' 
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Using  equations  1.3-18  and  1.3-19  it  is  seen  that 

P {t  <n5'<T-/p  } = 1-F(T  /n;p  )- 

r ( X ^ o f J.  o 1.3-21 

[l-F(T2/n;p^) ] 

or 

Pj.  |Tj^<n5'<T2/PQ  I i 1 - 2a  = Y 1.3-22 

Thus  the  probability  of  A being  in  the  'tails'  of  the  binominal 

distribution  is  less  than  2a  for  a fixed  p . 

_ _ o 

i.e.  , P (nA  <T. ) + P (nA>T_)  < 2a  1.3-23 

r — 1 r — 2 

The  two  contours  T^^,  T2  can  now  be  used  to  generate  confidence 
intervals  for  p with  confidence  coefficient  y. 

For  a particular  seunple  mean  K,  let  £ be  the 
largest  value  of  p such  that 

Tj^(Y/P)  5.  nA*  1.3-24 

and  let  p be  the  smallest  value  of  p such  that 


^2  (Y/P)  1 nA  1.3-25 

£ and  p depend  upon  the  saunple  mean  and  as  such  are  random 
variables  with  property  the  £ < p . 

£ is  called  the  lower  confidence  limit  and  p 
is  the  upper  confidence  limit.  The  interval  (£,p)  is  called 
the  100y%  confidence  interval.  This  interval  has  the  property 
that  if  p^  is  the  true  value  of  p,  then  the  probability  that 
the  interval  (£,p)  contains  p^  is  greater  than  or  equal  to  y. 

i.e.,  |e  1 Pq  1 P } i Y 1.3-26 

To  see  that  equation  1.3-26  is  valid,  consider  Figure  A-18. 
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Fiqure  A-18  Upper  and  Lower  Confidence  Limits 

A-44 
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T2  and  are  the  upper  and  lower  100y%  confidence  contours 
defined  by  equations  1.3-18  and  1.3-19.  Let  A be  the  average 
number  of  "successful"  hits  in  a random  sample  x^,  ...,  and 
assume  that  p^-  is  the  true  value  of  p.  Then  nA  is  the  number 
of  hits  above  1 in  this  sample  and  for  Tj^(Y»PQ)<nA<  T2(y»Pq) 
we  have 

^r  nA<T2(Y,P^)  /p^}  > Y 

Using  Figure  A-18  it  can  be  seen  that  (£,p)  contains  p^  if  and 
only  if 

Ti(y,Po)  1^1 
We  now  have 

Pr|Ti(Y»p)  1 nS”  1 T2(Y»P)/Pq  | 1 Y if  and  only  if 

Pr  1 Pq  - P } - 

Which  is  what  we  set  out  to  establish. 

Because  of  equation  1.3-12,  the  confidence 
interval  on  p can  be  related  to  one  on  p.  Let 


^ 


^ ln(p) 

If  p is  the  actual  mean  and  we  let 
o 
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then  we  have 

Pj.  {ii  i ~ 1 " { £ - Pq  - P I ^ Y 1.3-28 

This  gives  confidence  intervals  on  jj  with  confidence  coefficient 

Y. 


1,3.6  Determination  of  Confidence  Contours, 

(Y»P) » T2 (Y,P) 

It  can  be  seen  in  Subsection  1.3.4  that  from  a 
sample  ...,  of  receiver  noise,  computation  of  nA  gives 

the  number  of  hits  above  interval  1.  Using  A,  the  values  £ 
and  p or  jj  and  "y  can  be  calculated  for  a specific  confidence 
coefficient  y.  This  procedure  depends  heavily  upon  having  the 
functions  Tj^(y,p),  T2(Y/P)  for  0<p<l.  Calculation  of  T^  and 
T2  require  solving  equations  1.3-18  and  1.3-19.  Figures  A-13, 
A-16  and  A-17  were  obtained  by  a computer  program  evaluating 
equations  1.3-18  and  1.3-19  directly.  Figure  A-12  was  partially 
obtained  by  using  a normal  approximation  assuming  n is  large. 

As  n-x®  and  with  p fixed,  the  binominal  distribution  equation 
1.3-19  tends  to  the  normal  distribution  with  mean  np  and  variance 
np (1-p) . 

(x-np) ^ 

n— k _ 1 g 2np  ( 1— p ) 2^  3—29 

^2TTnp  (1-p)' 

Using  this  approximation  £ and  p can  be  determined  by  the 
(4) 

equation 


i.e., 


p(= 

1 

-T-.l.  A^a.  . 
nA  + ^ + A 

/(nA  + i)  (n-nA  + j 

E' 

n+  A*  a 

— 2 2 a 

V 

^ — 

1.3-30 

(4)  A.  Hald 
(Sixth 

, "Statistical  Theory  with 
printing  1965) , John  Wiley 
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Engineering  Applications", 
& Sons,  Inc. 
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where  is  the  solution  to  the  equation, 


a = 2 [1-0(A^)] 


a 


1.3-31 


The  normal  approximation  is  a "good"  approximation  if  np(l-p) 

M ) 

is  "large".  Error  bounds  can  be  obtained  for  this  approxi- 

mation and  in  the  calculations  of  ~ ,jj_  the  error  is  less  than 
3 X lO"^  if  nA  > 900. 


1 . 4 Equations  for  Computing  Weighting  Tables 

Weighting  tables  were  used  to  compensate  for  variations 
in  the  data  resulting  from  changes  in  range,  doppler,  and  angles 
off  antenna  boresight.  These  are  the  theoretically  known 
factors.  The  weighting  was  applied  to  the  output  of  the  FFT 
during  data  reduction.  Hardware  and  other  variabilities  were 
considered  separately  during  data  analysis. 


The  weighting  function  as  a multiplier  is  given  by 
the  expression: 


MULT(N,R,A) 


ANT(N,R,A)  ;AKEA(N,R,A)  . STTTC  (TT)  ! (A , R) 


' where:  N = doppler  cell  number 

R = range 
A = Altitude 

K = Dynamic  range  scaling 
ANT()  = Antenna  weighting 

AREAO  = Area  of  resolved  patch  on  ground 


(4)  A.  Hald,  "Statistical  Theory  with  Engineering  Applications", 
(Sixth  printing  1965),  John  Wiley  & Sons,  Inc. 

(5)  Harvard  Computation  Laboratory  (1955)  , "Tables  of  the 
Cumulative  Binomial  Probability  Distribution", 

Harvard  University,  introduction  pp.  XVIII 
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SINCO  = Sample/Hold  weighting 
0^0  = Grazing  angle  weighting 

This  expression  comes  directly  from  the  radar  range  equation, 
and  each  of  the  components  on  the  right  hand  side  of  the  equation 
is  described  as  follows: 

4 

a)  R Factor 

4 

The  computation  of  R is  a simple  function  of  the 
range  gate  number. 

b)  K Factor 

This  factor  is  chosen  to  optimize  dynamic  range  in 
the  succeeding  computations. 

c)  SING 

The  sin  x/x  frequency  aperture  generated  by  the  S/H 
was  obtained  by  simply  evaluating  the  sine  function 
at  the  appropriate  frequencies. 

d)  Mean  Radar  Cross  Section 

The  radar  cross  section  was  compensated  as  a 
function  of  grazing  angle  by  taking  data  from  the 
standard  curves.  This  gives  a fairly  good 
approximation  to  the  actual  function  and  reduces 
the  dynamic  range  required.  To  obtain  an  absolute 
measure  of  when  desired,  this  compensation  was 
reversed  in  the  data  analysis  procedure. 

e ) ANT 

The  antenna  pattern  (ANT)  equations  used  were  those 
in  Subsection  1.2.3  of  this  appendix.  A function 
of  the  crab  angle  estimated  from  the  data  during 
data  reduction  was  used  for  weighting.  Using 
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the  coordinate  system  of  Figure  A-19  the  following 
definitions  are  applicable: 

ANUM  - doppler  cell  number  starting  from  zero  doppler 
PRF  - sampling  rate  (14900  Hz) 

M - number  of  points  iii  FFT  (128) 

V - aircraft  velocity  (420  F/S) 

A - mean  wavelength  of  transmitter  output  (0.1  FT) 

G - antenna  gain  (lldB  on  axis) 

P - transmitter  power  (60dBm) 

RW  - range  dimension  of  equivalent  rectangular 

resolution  cell  on  sea 

a - angle  of  resolution  cell  off  aircraft  broadside 
in  horizontal  plane 

R - range  to  approximate  center  of  resolution  cell 

N - corrected  mean  value  of  receiver  noise 

The  coordinates  of  the  LOS  in  the  antenna  freime  are 
given  by: 

R = V cos  0 = C cos  'p  cos  p 
R = V sin  a cos  p 

R = J ^^(ANUM  - 0.5) 

„ . sin-1  PRF  (AHUM-0.5) 

2VM  cos  p 

tan  p = cos  a tan  8 = tan  ^ ^ ) 

oos  ex 

tan  0 = tan  ^ --  ) 

cos  a 

0 = tan  ^ (tan  a)  (cos  6) 
ip  = ip  ANT-tan'^  ( ) 
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V and  H define  plane  parallel  to  horizontal 
ground  plane 

I 

i 

s’  = proj  R on  vertical  plane  containing  iT 


Figure  A-19  Relative  Position  of  Aircraft  and 

Arbitrary  Doppler  Cell  on  the  Surface 
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f ) AREA 

The  area  of  each  range-doppler  cell  was  computed 
as  explained  below.  The  width  of  the  doppler  cell, 
i.e. , the  range  dimension  (x)  is  given  in  Figure  A-20. 


Figure  A-20  Vertical  Plan  Containing  the  LOS  to 
the  Leading  Edge  of  a Range  Gate 
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Then  the  length  of  a doppler  cell  is 

R = V cos  0 = V cos  ip  cos  p 
R = V sin  a cos  p 


^ _ 2R  _ 2V  cos  p sin  ot 

■ -J-  X 

2V 

6f^  = — cos  p cos  a d a 


da  = 


6f,  A 
a 


2V  cos  p cos  a 


6f,  AR_ 

_ j __  d (j 

~ 2V  cos  p cos  a 


R = 


cos  p 
PRF 


M 


A (PRF ) R 

D = M cosa  ' area  of  a doppler  cell  is 


A - 


A X PRF  X R X RW 
2 X V X N X cos  a 


1.5  Conditional  Probability  Maps  as  Diagnostic  Tools 

Conditional  probability  maps  were  orginally  intended 
to  characterize  the  neighborhoods  of  large  hits  and  were  so 
used  in  the  data  analysis.  These  maps  have  also  proved  to  be 
excellent  detectors  of  anomalous  behavior  in  runs.  This  property 
of  conditional  probability  maps  is  illustrated  below  using  run 
1603  (flight  16  run  3)  as  an  example. 

The  hit  count  versus  time  map  of  1603  (Figure  A-21) 
shows  a sudden  increase  in  hit  level  in  range  gate  2 at 
approximately  12.5  seconds.  This  unusually  high  level  con- 
tinues until  17.5  seconds. 
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Conditional  probability  maps  of  all  the  clutter  data 
in  flight  1603  are  shown  in  Figures  A-22,  A-23,  A-24,  A-25 
and  A-26.  That  the  phenomena  is  extended  in  time  is  readily 
seen  in  the  relative  range  = 0 map  (Figure  A-23)  and  the 
relative  doppler  = 0 map  (Figure  A-24) . From  the  time  collapsed 
array  (Figure  A-25)  it  is  clear  that  adjacent  doppler  cells  are 
involved.  From  the  normalization  array  (Figure  A-26)  it  can 
be  seen  that  the  third  range  gate  from  the  edge_of  the  maps 
contained  the  large  majority  of  hits  above  T^^.  This  is  of 
course  confirmed  by  the  hit  time  plot  (Figure  A-21) . 

Further  investigation  showed  that  erroneous  responses 
were  to  be  found  in  range  gate  2 (probably  due  to  computer 
problems  in  data  reduction) . Range  gate  2 was  deleted  from 
the  data  to  be  processed  and  the  resulting  conditional  prob- 
ability maps  (Figures  A-27,  A-28,  A-29  and  A-30)  show  the 
usual  (vertical  polarization)  structure  consisting  of  singular 
hits  at  the  0,0  point.  The  normalization  array  (Figure  A-31) 
shows  a biasing  of  the  large  hit  towards  one  side  of  the  map 
but  not  as  severe  as  in  Figure  A-26.  This  is  a direct  result 
of  deleting  range  gate  2 data. 

Conditional  probability  maps  have  proven  a valuable 
tool  in  locating  sometimes  subtle  processing  or  data  problems. 
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Figure  A-22  Conditional  Probability  with  Anomaly ,Time=0 
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Figure  A-23  Conditional  Probability  with  Anomaly, Range=0 
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Figure  A-25  Conditional  Probability  with  Anomaly , Time  Collapsed 
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Figure  A-26  Normalization  Array  with  Anomaly 
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Figure  A-28  Conditional  Probability  with  Correction , Range=0 
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Figure  A-29  Conditional  Probability  with  Correction , Doppler=0 
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APPENDIX  B 

OUTPUTS  - SPECIAL  ANALYSIS 


This  appendix  gathers  together  the  outputs  of  several 
special  analyses.  These  are  in  addition  to  the  standard  histo- 
gram and  other  outputs  which  were  done  as  a regular  part  of  the 
data  analysis.  Included  in  this  appendix  are  sections  and  sub- 
sections dealing  with  curve  matching/  histograms  of  the  mean, 
variation  of  statistics,  a special  case  of  a ship  in  the  field 
of  view,  behavior  in  the  vicinity  of  a large  hit  and  locally 
normalized  histogreuns  and  associated  plots. 

1.1  Curve  Matching  of  Histograms 

This  section  is  concerned  with  fitting  the  histogreuns  of 
the  sea  clutter  returns  to  some  mathematically  manageable  function. 
The  hope  is  to  find  an  expression  which  has  the  property  that 
the  difference  between  it  and  the  histogram  is  well  within  the 
Scunpling  errors  involved.  When  the  histogreun  data  was  first 
observed  there  were  attempts  to  characterize  the  data  as  coming 
from  a time  process  whose  mean  was  slowly  varying  and  whose 
first  Oder  distributions  were  Rayleigh  or  Weibull.  These  fits 
did  not  provide  the  best  match  since  they  did  not  track  very 
well  over  the  most  useful  part,  the  tails  of  the  observed  histo- 
grams. The  need  to  find  a curve /especially  one  which  fit  the 
tails, led  to  a match  by  cubic  splines.  Descriptions  of  the 
three  types  of  curve  matching  follow. 

1.1.1  Curve  Matching  Weibull  and  Rayleigh  Functions  to  Data 

There  are  many  theoretical  models  which  are  utilized 

as  statistical  models  of  clutter.  The  Weibull  model  which  has 
the  Rayleigh  model  as  a special  case  (N=l)  was  chosen  as  a 
possible  fit  for  the  clutter  data.  A least  squares  fit  of  the 
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Weibull  model  to  the  data  was  performed  and  graphical  results  of 
these  fits  are  given.  The  results  of  the  fits  show  the  Weibull 
model  to  be  a fair-to-poor  model  expecially  along  the  tails  of 
the  distributions.  For  the  Rayleigh  model  a least  squares  fit 
was  not  performed  but  the  mean  at  the  exponential  distribution 
was  made  equal  to  the  sample  mean  of  the  data.  For  this  model 
the  fit  was  worse  than  for  the  more  general  Weibull. 

The  results  of  the  fit  for  both  models  are  shown  only 
for  mean  adjusted  data  of  flight  605  with  vertical  polarization, 
and  flight  707  with  horizontal  polarization.  Figures  B-1  and 
B-2  are  plots  of  both  the  data  and  the  Weibull  fitted  for  runs 
605  and  707  respectively.  Figures  B-3  and  B-4  are  plots  of  the 
data  and  Weibull  model  with  axes  which  make  a Weibull  distribution 
a straight  line.  The  Weibull  distribution  is  given  by  the 
equation 

( 1 - EXP  - [— ] 0<x<“  n,  a>0 

F(x)  = { “ 

i 0 X < 0 

This  can  be  transformed  into  the  equation 


which  is  a straight  line  if  the  coordinate  axes  are 

ln[-  ln[-  (l-F(x)]]  versus  10  • 

The  least  squres  fit  of  a straight  line  to  the  data  was 
performed  in  this  coordinate  system  and  values  of  n and  x obtained 
are  given  in  these  graphs.  The  mean  may  be  related  to  a and  n 
by  the  expression 

X = p(l  + i ) 
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Figure  B-1  Data  Plot  and  Weibull  Fit  of  Log  Q From  Flight  605 
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Figure  B-2  Data  Plot  and  Weibull  Fit  of  Log  Q From  Flight  707 
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Figure 

B-4  Data  Plot  and  Weibull  Fit  of  Flight 
707  in  Weibull  Coordinates 
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A measure  of  the  goodness  of  fit  is  given  by  graphs  B-5, 
B-6,  B-7  and  B-8. 

For  run  605  it  can  be  seen  that  from  the  10  ^ to  10  ^ 
points  the  model  can  differ  from  these  points  by  as  much  as 
0.25dB.  For  707  this  difference  is  high,  about  0.75dB.  However, 

-3 

if  the  ratio  of  the  mean  to  the  10  point  is  considered,  the 

error  can* be  2dB.  As  can  be  seen,  the  data  from  the  vertical 

•» 

polarization,  flight  605,  is  fitted  better  by  the  Weibull  model 
than  flight  707  with  horizontal  polarization.  In  either  case, 
the  Weibull  is  not  an  especially  good  fit  when  compared  to  the 
cubic  spline  fit  of  the  next  subsection. 

The  results  for  the  Rayleigh  fit  are  worse  than  those 
for  the  Weibull.  One  need  only  look  at  Figures  B-9  and  B-10 
to  see  the  poorer  fit.  The  tails  of  the  exponential  distribution 
are  much  lower  than  the  tails  of  the  data.  The  errors  in  the 
match  are  given  by  Figures  B-11  through  B-14.  Errors  with  respect 
to  the  matched  mean  are  1.8dB  at  the  10  ^ point  for  run  605  and 
3.8dB  for  run  707. 

If,  rather  than  attempting  a Weibull  least  squares  fit, 

we  try  only  for  a match  at  the  mean  and  a close  fit  in  the 
—3  *"6 

Q = 10  through  10  region  (but  still  use  the  Weibull  plots) 

the  problem  becomes  easier.  The  reason  for  doing  this  would  be 

related  to  the  usual  radar  CFAR  use  where  the  mean  (or  a point 

reasonably  close)  may  be  readily  measured  and  false  alarm 

thresholds  can  then  be  selected  based  on  an  extrapolation  to 
— 3 —6 

the  10  or  10  point.  Pragmatism  rather  than  rigor  is  the 
guide  and  a fit  in  these  terms  to  better  than  l/2dB  is  quite 
acceptable.  For  this  case,  the  fit  for  run  605  is  already 
acceptable  and  run  707,  an  extreme  case,  almost  satisfies  the 
l/2dB  criterion.  Figure  B-15  is  simply  Figure  B-4  with  a new 
Weibull  fit  from  the  mean  through  +15dB  (corresponding  to 
Log  Q = -6.0).  This  heuristic  fit  from  the  mean  through  the 
10  ^ point  is  within  0.75dB. 
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Figure  B-6  Error  of  Weibull  Fit  to  Q(x)  Data  From  Flight  605 
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Figure  B-7  Error  of  Weibull  Fit  to  Log  Q Data  From  Flight  707 
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Figure  B-8  Error  of  Weibull  Fit  to  Q(x)  Data  From  Flight  707 
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Figure  B-9  Data  Plot  and  Rayleigh  Fit  of  Log  Q From  Flight  605 
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• • B-IO  Data  Plot  and  Rayleigh  Fit  of  Log  Q From  Flight  707 
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Figure  B-11  Error  of  Rayleigh  Fit  to  Log  Q Data  From  Flight  605 
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Figure  B-12  Error  of  Rayleigh  Fit  to  Q(x)  Data  From  Flight  605 
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Figure  B-14  Error  of  Rayleigh  Fit  to  Q(\)  Data  From  Flight  707 
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1.1.2  Cubic  Spline 

The  cubic  spline  fit  was  made  to  the  histogram  of  the 
data  "x"  expressed  as  a function  H which  is  related  to  the 
histogram  through  the  observed  cumulative  distribution.  The 
function  H is  defined  as  h(x)  = QM  where  Q(x)  = l-F(x) 

and  F(X)  is  the  cumulative  distribution  function  of  the  observed 
variable  x. 

The  function  log^ - Q has  the  advantage  that  it  explicitly 
-k 

gives  the  10  points  k = 1,2,3,...,  and  also  gives  a better 
description  of  the  tails.  The  values  are  identified  by  the 
equation 


P (x  > a,  )—  . k — 0,1,2,.... 

r K X 

Using  this  definition  of  we  can  easily  see  that 

H(aj^)  - logj^QQ(Aj^)  = log^Q[10“^]  = -k.  k = 0,1,2,.... 

Thus  a value  of  -5  for  logj^QQ(a)  means  that  "a"  is  that  number 

such  that  the  probability  of  obtaining  any  value  greater  than 

it  is  — ^ . The  usefulness  in  having  such  points  directlv 
10^ 

available  is  in  easily  observing  important  threshold  levels 
especially  those  needed  in  determining  false  alarm  rates. 

a)  Curve  Matching  to  Cubic  Spline 

The  spline  function  is  defined  in  terms  of  a 

set  of  points  called  knots.  If  the  m numbers  tk^^  < tk2<.*.<tkj^ 

are  a set  of  knots  then  the  spline  S(t)  is  a polynominal  of 

degree  n between  the  knots  with  the  condition  that  it  have  n-1 

continuous  derivatives  at  the  knots.  Specifically  for  n=3  the 

spline  S(t)  is  called  a cubic  spline  and  it  is  a piecewise 

cubic  polynominal  over  the  interval  (tkj^,tkjj^)  which  has  the 

property  that  it  is  continuous  and  has  continuous  first  and 

second  derivatives.  For  knots  tl,  < tk-  < ...  < tk  the  cubic 

i.  z m 

spline  S(t)  is  usually  expressed  as 
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S(t)  - c^_3  + c._3  D + 

with  D = t - tk^  and  tk^  f.  ^ ^ ^^i+1" 

Complete  descriptions  then  of  a cubic  spline 
requires  knowing  the  knot  values  and  for  each  pair  of  consecutive 
knots  the  four  coefficients  in  the  cubic  equation.  The  description 
of  a cubic  spline  approximation  to  the  data  given  in  this  report 
does  not  follow  the  usual  convention.  It  was  decided  to  use 
the  format  that  if  S(t)  is  cubuc  spline  approximation  to  H(x) 
where  t = log^^pCx/lO)  with  four  knots  tk^  < tk2  < tk^  < tk^ . 


Then  for  values  tk^  1 t 
can  be  written  in  the  form 


i = 1,2,3,...  so  that  S 


S(t)  = J . A._3  t . A._^ 

Where  t = 10  log^^^  (x/10)  . 

The  above  format  of  S(t)  avoids  the  intermediate  variable  D 
and, for  a complete  description  of  a cubic  spline  with  four  knots, 
we  need  only  a 3 X 4 matrix  A defining  the  coefficient  of  the 


and  a 

matrix 

tk  of 

knot 

1,1 

^1,2 

^1,3 

^,4 

2,1 

^2,2 

^2,3 

^2,4 

3,1 

^3,2 

^3,3 

^3,4 

, tk3 

, tk^) . 

The 

matrix 

the  graph  of  the  cubic  spline  approximation  for  each  data  set 
fitted. 

For  fixed  knots  the  cubic  spline  approximation, 
S(t),  to  the  data  H(Xj^)  i = 1,2,...,  N is  made  through  a least 
square  fit.  The  quantity,  E,  which  is  minimized  is 
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E = 


N 


[H(X^)  - S(t^)] 

Li=l 


1 

2 


where 


t.  = 
X 

10  • 

logio(Xi/10) 

i = 1,2, . . . ,N 

W.  = 
1 

w.  = 
1 

(^i+1 

ti_i)/tN  - 

i ^ 2,3,* 

«N  = 

- 

The  key  to  the  success  of  a spline  approximation 
is  the  location  of  the  knot  values.  A variable  knot  computer 
progrcun  was  used  to  determine  in  an  economical  fashion  the 
placement  of  the  knots  which  gives  the  minimum  least  squares 
error.  The  computer  progreun  is  a standard  routine  in  the 
International  Mathematical  and  Statistics  Libraries  (IMSL) 
called  ICS  VKU. 


b)  Svunmary  of  Cubic  Spline  Match 

Cubic  splines  of  various  degrees  were  matched 
to  the  data  and  it  was  found  that  the  use  of  four  knots  was 
sufficient  to  approximate  down  to  the  sample  noise  of  the 
data.  In  this  section  we  give  only  results  on  the  fit  by  cubic 
splines  with  only  four  knots. 

I 

The  chart  on  Table  B-1  gives  the  flights  and  runs  | 

that  were  fitted  by  splines.  There  was  an  attempt  to  fit  a ' 

good  cross  section  of  the  available  data.  Flight  707  and  1706  « 

were  the  only  flights  with  horizontal  polarization  fitted  and 
flight  803  was  a low  grazing  angle  run.  The  data  from  flight 
605  were  processed  two  different  ways,  one  as  A type  and  once 
as  N type.  The  process  involved  in  forming  the  A type  and  N type 
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TABLE  B-1 

FLIGHT  INFORMATION 
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is  described  in  detail  in  Volume  II.  The  difference  between 
them  is  that  the  N type  was  an  attemp  to  remove  the  variations 
in  mean  clutter  levels  that  were  noticed  during  the  runs.  Cubic 
splines  approximations  were  made  to  the  two  different  processes 
of  the  data. 

For  each  flight  given  in  Table  B-1  there  corres- 
ponds four  graphs.  The  four  graphs  for  each  flight  are  titled 
as  follows: 

1 . Log  Q 

2.  Cubic  Spline  Approximation  to  Log  Q 

3.  Error  of  Fit  of  S(t)  to  Log  Q 

4.  Error  of  Fit  of  10®^^^  to  Q 


The  first  graph,  ’!og  Q,  is  the  actual  data  which  is  being  matched. 
The  second  graph,  Cubic  Spline  Approximation  to  Log  Q,  is  a 
graph  of  the  cubic  spline  approximation.  The  final  knot  locations 
are  identified  by  x's  marked  on  this  curve  and  matrices  are  given 
identifying  the  cubic  spline  that  is  plotted.  In  a)  above 
a description  is  given  on  how  to  construct  the  cubic  spline  from 
the  matrix  A and  tk  which  are  given  below  each  cubic  spline  curve. 

The  third  graph.  Error  of  Fit  of  S(t)  to  Log  Q, 
is  a plot  of  the  difference  D between  Log  Q(x),  the  data  set, 
and  S(t),  t = 10  • log^^O  Yq  cubic  fit,  i.e.  , D(t)  = 

Log  Q (10^^  - S(t).  Along  with  the  plot  the  least  squares 

error  of  the  cubic  spline  approximation  is  given.  This  graph 
gives  us  a measure  of  how  well  the  cubic  spline  fits  Log  Q. 

One  can  see  that  in  most  cases  the  cubic  spline  fits  out  to  the 

-4  -4 

10  point  with  an  error  of  at  most  .02.  Beyond  the  10  point 

the  error  curve  indicates  that  the  spline  is  a fit  down  to  what 

appears  as  Scimpling  noise  due  to  the  smaller  number  of  samples 

at  these  levels. 
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The  fourth  graph  of  a set.  Error  of  Fit  of 
loS(t)  to  Q is  a measure  of  the  fit  of  the  spline  surve  in  terms 
of  Q.  The  quantity  graphed,  E (b) , is  given  by  the  equation 

E(t)  = QdO-^  nr  ) - S(t) 

where  S(t)  is  the  cubic  spline  fitted  to  Log  Q.  The  horizontal 
axis  of  all  graphs  given  are  in  terms  of  t in  dB)  rather 

than  in  terms  of  x directly.  The  relation  between  t and  the 
observed  variable  x is 

t = 10  [log^^Q  (x/10)  ] 

c)  Graphical  Results  of  Cubic  Spline  Approximation 
The  succeeding  figures  are  the  graphical  results 
of  the  spline  curve  fit.  To  aid  in  finding  a particular  graph 
the  cross-reference  of  Table  B-2  may  be  used.  This  charts 
tabulates  the  figure  numbers  for  the  type  of  graph  and  data  set 
desired. 


TABLE  B-2 

FIGURE  NUMBER  LOCATOR 


605-A 

403-N 

605-N 

707-N 

803-N 

1104-N 

1602-N 

1706-N 

LOG  Q 

B-16 

B-20 

B-24 

B-28 

B-32 

B-36 

B-40 

B-44 

Cubic 
Spline 
Approx  to 
LOG  Q 

B-17 

B-21 

B-25 

B-29 

B-33 

B-37 

B-41 

B-45 

Error  of 
Fit  S(t) 
to  Log  Q(x) 

B-18 

B-22 

B-26 

B-30 

B-34 

B-38 

B-42 

B-46 

Error  of 
Fit  to 

10="=’ 
to  Q(x) 

B-19 

B-23 

B-27 

B-31 

B-35 

B-39 

B-43 

B-47 

In  summary,  it  may  be  simply  stated  that  there  is 

a Maximum  error  between  the  cubic  spline  fit  and  the  data  of  less 
than  O.ldB.  ^ 
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Fig.  B-17  Cubic  Spline 

Approximation  to 
LOG  Q's  Run  605A 
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Figure  B-21  Cubic  Spline 
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Fig.  B-25  Cubic  Spline 

Approximation  to 
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Figure  B-26  Error  of  Fit  of  S(t) 
to  Log  Q,  Run  605 
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1.2  Histograms  of  the  Mean 

The  mean  was  determined  by  averaging  over  all  range 
gates  and  doppler  cells  and  10  freunes  in  time  (i.e.,  5120  data  * 
points) . There  are  approximately  3000  such  means  on  a typical 
5 minute  data  run.  Reduced  number  of  range  gates  or  time 
duration  in  some  runs  reduce  the  number  of  seunples  per  run. 

These  Scunple  means  result  in  data  streeun  of  approximately  3000 
scunples  with  a seunple  rate  of  11.1  means/  second. 

A simple  histogreun  was  found  from  these  mean  values  over 
a run.  These  histograms  in  pdf  form  are  displayed  in  the 
succeeding  pages  for  selected  cases.  The  temporal  correlation 
of  the  mean  (from  autocorrelation  analysis)  is  high  so  that  the 
number  of  independent  mean  Scunples  is  low  even  for  a 5 minute 
run.  Some  of  the  irregular  structure  in  these  histograms  comes 
from  this  limited  data  base.  Much  longer  data  runs  would  be 
required  to  obtain  data  bases  sufficient  for  more  extensive 
mean  analysis. 
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1 . 3 Temporal  Spectra  and  Autocorrelation  Function  of  the  Mean 

The  data  set  for  spectra  and  correlation  analysis  is 
the  same  as  for  histograms  of  the  mean  of  the  precesding  sub- 
section. A 1024  sample  selection  of  the  data  was  Fourier  trans- 
formed to  obtain  a spectrum.  Real  and  imaginary  parts  were 
squared  and  summed  to  obtain  the  PSD.  The  PSD  was  then  Fourier 
transformed  to  obtain  the  autocorrelation  function  (ACF) . This 
transform  was  aperiodic  and  compensated  to  obtain  the  ACF.  The 
aperiodic  compensation  and  higher  sensitivity  produced  increased 
noise  on  the  tails  of  the  ACF. 

Again  here  the  limited  number  of  independent  mean  samples 
creates  some  of  the  long  term  oscillation  in  the  ACF.  This 
assessment  has  been  born  out  by  the  simulation  results  (see 
Section  10  of  Volume  II  and  Appendix  F of  this  volume) . 
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1. 4 Variability  of  Selected  Statistics  j 

Selected  statistics  of  mean  normalized  histogreuns  were  I 

analyzed  for  variability. 

-1  -2  -3  -4 

The  variance  of  the  median/  10  / 10  , 10  and  10 

points  were  computed.  Both  the  mean  value  and  the  variance, 
organized  by  flight  and  run,  are  listed  in  Tables  B-3  through  B-9. 

Before  using  this  data  in  further  analysis  one  should  j 

check  hit  maps,  hits  vs.  time,  the  distribution  function  and 
scenarios  to  ensure  that  only  valid  and  comparable  data  is 
being  analyzed.  For  example: 

1.  Horizontal  polarization  runs  7 and  17  show 
considerably  more  variability  tham  most  vertical 
polarization  runs. 

-4 

2.  The  very  high  variability  of  the  10  point  versus 
its  value  on  812  is  worth  noting.  812  was  the  run 
containing  the  light  ship  and  as  such  describes  its 
characteriztics  rather  than  that  of  the  sea. 

3.  Two  run  403 's  are  shown,  403  (-1)  has  range  gate  1 

I 

deleted.  Range  gate  1 is  included  in  403 (All). 

Run  403  had  data  reduction  errors  in  range  gate  1.  i 
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TABLE  B-3 

VARIABILITY  OF  SELECTED  STATISTICS 
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TABLE  B-4 

VARIABILITY  OF  SELECTED  STATISTICS  (cont'd) 
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TABLE  B-5 

VARIABILITY  OF  SELECTED  STATISTICS  (Cont'd) 
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VARIABILITY  OF  SELECTED  STATISTICS  (Cont'd) 
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TABLE  B-7 

VARIABILITY  OF  SELECTED  STATISTICS  (Cont'd) 
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TABLE  B-9 

VARIABILITY  OF  SELECTED  STATISTICS  (Cont'd) 
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1 , 5 Hit  Map  of  Run  812  - Nantucket  Light  Ship  Target  Validation 

Flight  812  was  a special  test  in  which  a known  target  was 
present  in  the  sea  clutter  data.  The  target  was  a light  ship  off 
the  eastern  coast  of  Massachusetts  positioned  off  Nantucket  Island. 

A special  Hit  Map  computer  run  was  executed  to  validate 
the  synthetic  aperture  mapping  process.  The  threshold  was  set 
at  500  — well  above  the  clutter  level  so  that  target  returns 
from  the  ship  could  be  easily  distinguished.  The  flight  was  run 
at  500  feet  altitude. 

Hits  exceeding  threshold  for  the  target  were  dumped  from 
a computer  run.  The  results  are  plotted  in  Figure  B-48.  Range 
gate  number,  doppler  filter  number  and  FFT  time  frame  for  hits 
exceeding  threshold  were  plotted  to  show  the  range-doppler  time 
history.  The  target,  initially  observed  in  range  gate  5 and  at 
the  highest  doppler  filter  (51),  passed  through  four  range  gates, 
and  all  32  doppler  filters  at  essentially  a linear  rate  in  about 
770  FFT  frames,  or  7.04  seconds. 

If  the  doppler  rate  is  assumed  linear,  velocity  of  the 
aircraft  with  respect  to  the  target  may  be  estimated.  This 
information  is  required  for  motion  compensation. 


Doppler  frequency  is  related  to  velocity  and  look  angle 


by 

^d  = ® 

where  A is  radar  wavelength  (0.1  feet).  Flight  test  geometry 
is  indicated  in  Figure  B-49. 

Doppler  rate  is  computed  by  differentiation  of  doppler 
frequency  with  respect  to  time: 


dv  Q 

dt  ® 


V sin  0 


de 

dt 
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If  aircraft  velocity  is  constant  the  first  term  is  small 
and  will  be  neglected.  From  the  geometry  the  velocity  is  related 
to  angular  rate  (d6/dt)  by 


Solving  the  last  two  equations  for  velocity  gives  a simple  result 
- XR(df  ,/dt) 

2 sin  0 

From  range-doppler  history  the  doppler  rate  is  estimated 
from  number  of  doppler  filters  (32) , doppler  bandwidth  (14000/ 

128  = 109.375  Hz)  and  observation  time  noted  earlier  (7.04  seconds). 


Figure  B-49  Flight  Test  Geometry  (Plan  View) 

B-112 

UNCLASSIFIED 


UNCLASSIFIED 


df 


dt 


d _ 32  X 109.375 


7.04 


= -497.16  Hz/seconds. 


If  nominal  range  (R  = 4850  ft)  and  angle  (©  = 90®)  are 
assumed  the  velocity  for  the  X-Band  radar  ( A = 0.1  ft)  is 
computed  as 


V = 


-0.1  X 4850  X (-497.16)  _ 


2 sin  (90®) 


= 374  ft/sec. 


The  velocity  estimate  was  used  to  generate  the  hit  map 
printouts  in  Figure  B-50.  In  the  upper  map  hits  exceeding  the 
threshold  are  plotted  in  each  range  gate  vs.  down-range  distance. 
The  lower  map  shows  hits  mapped  into  a cross-range  vs.  down- 
range  coordinate  system  quantized  to  100  foot  cells.  The  range 
gate  sum  (when  totaled)  add  up  to  1215  hits  exceeding  threshold 
during  observation  of  the  Lightship  Target. 

1.5.1  Lightship  Conditional  Probability  Maps 

A complete  set  of  conditional  probability  maps  was 
generated  for  Flight  8 Run  12.  Contour  plots  exist  for  all  cuts 
in  time,  range  and  doppler  through  the  conditional  probability 
cube.  A time  collapsed  array  and  normalization  array  were  plotted 
as  well.  These  plots  validate  the  conditional  probability 
mapping  process  against  the  range-doppler-time  history  plotted  in 
Figure  B-48. 

The  conditional  probability  time  cuts  in  Figures  B-51 
through  B-58  show  the  Lightship  Target  is  extended  from  minus 
seven  to  zero  FFT  time  frames  (one  FFT  frame  = 0.0091428571  sec). 
The  range-doppler  history  indicates  the  target  is  well  extended 
over  several  FFT  time  frames  in  each  range-doppler  cell.  The 
conditional  probability  map  cut  at  relative  time  = 0 . 

(Figure  B-58)  shows  doppler  extent  is  a little  greater  than 
range  extent  for  the  Lightship  Target.  The  result  is  in 
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Figure  B-50  Hit  Maps  for  Flight  0812  Nantucket  Lightship  Target 
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agreement  with  range-doppler  history;  in  particular,  hits  exceed 
threshold  in  either  one  o"  two  range  gates  and  one  to  three 
doppler  cells  at  any  po.  in  time. 

Figures  B-59  to  B-62  display  conditional  probability 
contour  plots  for  relative  range  -1,  0,  +1,  +2.  All  other  range 
cuts  were  void  of  data  and  have  therefore  been  deleted.  Very 
few  points  are  indicated  outside  of  the  relative  range  = 0 cut; 
in  fact,  the  number  of  hits  was  not  enough  to  cause  the  first 
contour  level  to  be  exceeded.  Cuts  in  doppler  containing 
information  are  shown  in  Figures  B-63  through  B-65.  Cuts  at 
relative  dopplers  farther  removed  from  the  +1  points  from 
relative  doppler  equal  zero  had  very  few  hits.  Only  the  three 
cuts  shown  had  tiny  probability  contours  (^0.05)  which  could  be 
plotted. 

The  time  collapsed  array  in  Figure  B-66  shows  the  target 
is  extended  more  in  doppler  than  range.  The  range  extent  is 
never  greater  than  two  range  gates.  Total  doppler  extent  is 
about  three  cells. 

1.5.2  Lightship  Histograms 

Absolute  statistics  data  are  plotted  in  Figures  6'6? 
and  B-68.  The  first  plot  shows  absolute  amplitude  in  decibels 
plotted  against  range  gate  number.  Results  indicate  saturation 
in  range  gates  2-5,  in  agreement  with  the  range-dopplor  time 
history  in  Figure  B-48.  Saturation  occured  during  an  extended 
time  period  (4  histogram  frames  = 2400  FFT  frames)  in  Figure 
B-48.  The  cimplitude  peak  in  histogram  record  frame  five  correlates 
well  with  time  the  lightship  was  known  to  be  present.  Frames 
4 and  5 show  the  peak  which  corresponds  to  2400  to  3000  FFT 
frames  in  Figure  B-48. 
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Target  Conditional 
Figure  B-60  Probability  Map 
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1 . b Characteristics  in  the  Neighborhood  of  Large  Returns 

• Conditional  probability  maps  presented  elsewhere  in 

this  report  give  an  overall  picture  of  large  hits  and  present 
information  on  returns  in  their  immediate  space-time  vicinity. 

This  section  gives  another  view  of  the  same  phenomena. 

1.6.1  History  of  a Typical  Large  Sea  Echo 

The  nature  of  a large  echo  when  horizontal 

polarization  is  used  is  shown  in  Figure  B-69.  A bar  graph  of 

power  versus  time  has  been  plotted  for  each  doppler  filter.  Bars 

were  used  since  the  output  was  a series  of  values  occurring  as 

every  FFT  was  completed.  Only  those  values  above  the  threshold 
2 

(8m  ) in  range  gate  3 are  plotted.  The  patch  of  sea  started 
ringing  filter  34  first.  Then,  as  the  aircraft  flew  past  the 
patch,  the  doppler  shift  was  reduced  as  a consequence,  causing 
the  echo  to  pass  through  several  filters  in  succession.  (In 
the  diagram  an  increasing  doppler  filter  number  corresponds  to 
a decreasing  doppler  frequency.) 

It  can  be  seen  that  the  output  of  each  filter  is 
fluctuating,  rather  than  constant,  and  that  during  the  time 
shown  the  probability  of  crossing  the  threshold  is  quite  high 
compared  with  the  rest  of  the  hit  map.  In  fact  there  are  no 
returns  above  the  threshold  which  are  not  on  the  main  diagonal. 
This  suggests  that  this  series  of  large  echoes  come  from  a 
relatively  stationary  part  of  the  sea  where  the  mean  cross- 
section  was  temporarily  high  and  lasted  for  approximately  1.5  sec. 
The  patch  extended  over  three  doppler  cells  at  its  maximum, 
implying  an  extent  of  one  or  two  hundred  feet. 

More  than  one  large  response  may  be  seen  in 
some  of  the  doppler  filters  (e.g.  39).  This  causes  the  conditional 
probability  map  to  extend  in  time,  causing  the  pattern  seen  in 
Figure  B-70. 
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Figure  B-69  History  of  a Typical  Large  Sea  Echo 
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!• 7 Unconditional  vs.  Conditional  Probability  Analysis 

Detailed  computer  analysis  performed  on  conditional 
probability  cubes  indicates  conditional  probability  of  getting 
a hit  in  the  vicinity  of  a large  hit  is  greater  than  unconditional 
probability.  This  property  is  illustrated  in  the  bar  graph  in 
Figure  B-71.  This  result  is  not  unexpected  - large  clutter 
returns  generally  occur  in  clumps  and  large  returns  have  associated 
relatively  large  returns  in  the  immediate  space-time  vicinity. 

Unconditional  probability  was  computed  from  the  ratio 
of  total  hits  exceeding  threshold  (from  the  Hit  Map  printout)  to 
total  number  of  hits  possible  for  the  run  (Histo  records  X 600  FFT 
fraunes/record  X number  of  doppler  filters  X number  of  range 
gates)  . 

Conditional  probability  was  computed  from  the  ratio  of 
total  hits  to  total  trials  within  the  conditional  probability 
cube.  The  central  cell  in  the  cube  was  excluded  from  the 
calculation. 

Conditional  probability  exceeded  unconditional  probability- 
in  all  cases  except  0801.  The  simulation  results  are  included 
in  the  results  at  the  bottom  of  the  Figure  B-71. 
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1.8  Locally  Normalized  Histogreuns 

The  histograms  in  this  section  are  designated  TOTAL  N. 

They  differ  from  the  histograms  in  Volume  IV  in  that  they  are 
locally  normalized.  Each  clutter  histogram  transferred  from  data 
reduction  is  normalized  by  its  mean  before  combining  to  obtain 
a total  histogr5un. 

N Histogreuns  were  obtained  for  both  typical  runs  and  for 
the  runs  with  extreme  tails  in  the  A histogreuns. 

These  histogreuns  are  defined  mathematically  as  TOTAL 

r.  . H.  . where  H . . is  H. . normalized  by  the  mean  of  H. ..  The 
13  i3  1-3  1-3  13 

histogreuns  are  displayed  in  each  of  the  standard  output  forms. 

For  each  output  form  the  histogreuns  are  arranged  sequentially 
by  flight  and  run. 

1.8.1  Histograms  TOTAL  N LOG  Q 

All  valid  clutter  data  within  each  run  is  included 
in  the  histogreun. 

The  N indicates  that  each  of  the  histogreuns  H^^ 
was  normalized  by  its  mean  before  combining  to  form  the  total 
histogreun  for  the  run.  The  vertical  axis  is  the  logarithm  base 
10  of  1 minus  the  cumulative  probability,  i.e.,  LOGj^q  (1-P  ) . 

Various  points  on  the  tail  of  the  distribution  are  clearly  read 
from  this  type  of  plot.  The  horizontal  axis  is  the  clutter  power 
per  cell  in  dB  referenced  to  the  data  mean. 

This  type  of  outputs  is  described  in  more  detail 
in  Section  9 (Volume  II) . Procedures  utilized  in  developing 
the  histogreun  from  the  raw  data  are  described  in  Section  8 of 
Volume  II. 
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1.8.2 


All  valid  clutter  data  for  each  run  is  included. 


The  suffix  (N)  after  TOTAL  indicates  each  histo- 
greun  was  normalized  by  its  mean  before  combining  into  the  total 
histogram  for  the  run.  The  vertical  axis  of  Weibull  plots  is 

In  (In ( ) where  P(x)  is  the  cumulative  probability  density 
function.  The  horizontal  axis  is  the  clutter  power  per  cell  in 
dB  referenced  to  the  data  mean. 

This  type  of  output  is  described  in  more  detail 
in  Section  9 (Volume  II) . Procedures  utilized  in  developing  the 
histogram  from  the  raw  data  are  described  in  Section  8 of  Volume  II, 
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1.8.3  Histograms  TOTAL  N Q = 1-DIST 

All  valid  clutter  dta  for  each  run  is  included. 

The  suffix  (N)  after  TOTAL  indicates  each  histo- 
gram was  normalized  by  its  mean  before  combining  into  the  total 
histogram  for  the  run.  The  vertical  axis  is  Q or  one  minus  the 
cumulative  distribution.  The  horizontal  axis  is  the  clutter 
power  per  cell  in  dB  referenced  to  the  data  mean. 

This  type  of  output  is  described  in  more  detail 
in  Section  9 (Volume  II) . Procedures  utilized  in  developing  the 
histogreun  from  the  raw  data  are  described  in  Section  8 of  Volume  II. 
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1.8.4  Histograms  TOTAL  N PDF 

All  valid  clutter  data  for  each  run  is  included. 

The  suffix  (N)  after  TOTAL  indicates  each  histo- 
greun  was  normalized  by  its  mean  before  combining  into  the  total 
histogram  for  the  run.  The  vertical  axis  is  the  raw  histogreun 
data  or  probability  density.  The  horizontal  axis  is  the  clutter 
power  per  cell  in  dB  referenced  to  the  data  mean. 

This  type  of  output  is  described  in  more  detail 
in  Section  9 (Volume  II) . Procedures  utilized  in  developing 
the  histogram  from  the  raw  data  are  described  in  Section  8 Volume  II. 
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1.8.5  Histoqreuns  TOTAL  N PDF  TAIL 

All  valid  clutter  data  for  each  run  is  included. 

The  suffix  (N)  after  TOTAL  indicates  each  histo- 
graun  was  normalized  by  its  mean  before  combining  into  the  total 
histogram  for  the  run.  The  vertical  axis  is  probability  density. 
Only  the  tail  is  displayed  to  give  dynamic  range.  The  horizontal 
axis  is  the  clutter  power  per  cell  in  dB  referenced  to  the  data 
mean. 

This  type  of  output  is  described  in  more  detail 
in  Section  9 (Volume  II) . Procedures  utilized  in  developing  the 
histogreun  from  the  raw  data  are  described  in  Section  8 Volume  II. 
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APPENDIX  C 

ENVIRONMENTAL  CONDITIONS 

Determination  of  wind  and  sea  conditions  was  done  as 
part  of  the  flight  test  program  and  was  a General  Duneunics 
responsibility.  The  following  material  was  excerpted  from  their 
flight  test  report  and  is  included  here  for  completeness  and 
for  the  convenience  of  the  reader. 

Sea  and  wind  conditions  and  the  resultant  sea  state 
were  recorded  for  each  flight  as  summarized  in  Table  C-1. 

Surface  contacts  noted  were  also  recorded  (Table  C-2) . 

More  complete  weather  data  over  a period  of  hours 
before  most  flights  was  also  recorded  and  is  repeated  from 
the  General  Dyneimics  flight  test  report  as  Table  C-3  through 
C-17. 
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TABLE  C-1 

SEA  CONDITIONS  FOR  TAGSEA  CLUTTER  FLIGHTS 
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TABLE  C-2 

SUIVIMARY  OF  SEA  SURFACE  CONTACTS 


FLIGHT 

NO.* 

Date 

Runs  With  Pictures  of 
Sea  Surface  Objects 

Frame 

Nos. 

Location  of 
Sea  Surface  Objects 

Approx.  Time 
T = Start  of  Run** 

2.2  kft  cross  wind 

25 

None  within  radar 
range  gate  coverage. 

3E 

2/5/76 

3.3  kft  up  wind 

44 

None  within  radar 
range  gate  coverage. 

3.3  kft  down  wind 

23 

None  within  radar 
range  gate  coverage. 

1 

Ship  approx  in  range 
gates  10,  11,  12,  & 13 

T + 7 sec. 

6E 

2/12/76 

2.2  kft  down  wind 

5 

None  within  radar 
range  gate  coverage 

15 

Ship  approx  in  range 
gates  14,  15,  & 16. 

T + 1 min;  45  sec. 

2.2  kft  up  wind 

59 

Ship  appprox  In  range 
gates  14,  15,  & 16. 

T 6 min:  53  sec. 

8 

None  within  radar 

7E 

. 

2/12/76 

2.2  kft  cross  wind 

9 

range  gate  coverage 

None  within  radar 
range  gate  coverage 

IIW 

3/10/76 

64 

Ship  approx,  in  range 
gates  3 and  4. 

T 7 min;  28  sec. 

65 

Ship  approx.'  in  range 
gates  1 and  2. 

T + 7 min:  35  sec. 

12W 

3/12/76 

2.2  kft  down  wind 

54 

None  within  radar 
range  gate  coverage 

* E • East  Coast 
W = West  Coast 

ONLY  FLIGHTS  CONTAINING  PICTURES 
OF  SEA  SURFACE  OBJECTS  ARE  LISTED. 

• • 

EXACT  IRIG  TIME  CAN  BE  DETERMINED  BY 
CORRELATING  FRAME  NUMBERS  FROM  EACH 
RUN  TO  CAMERA  FRAME  PULSES  (TRACK  3 
OF  MAGNETIC  TAPE  RECORDER). 
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APPENDIX  D 

DATA  ANALYSIS  SOFTWARE 


This  appendix  describes  the  software  used  in  TAGSEA  data 
analysis.  Functionally  the  tasks  performed  by  the  software  are 
described  in  Section  8 (Vol.  II). 

Software  structure  for  histograms  and  hit  processing  are 
documented  here.  Software  for  average  analysis  is  included  with 
that  for  histograms  since  this  was  the  way  the  software  was 
structured  as  some  commonality  exists.  Histogram  software  is 
divided  into  two  parts,  initial  and  post  processing,  again  based 
on  the  way  the  software  was  written. 

1 . 1 Histogram  and  Mean  Data  Analysis  Software 

Several  programs  are  used  to  perform  histogram  and  mean 
analysis.  The  interrelationship  of  the  program  is  shown  in 
Figure  D-1. 

A description  of  each  of  the  programs  follows. 

1.1.1  HISTOGRAM  PROGRAM  (TAGHIST) 

The  computer  progreim  TAGHIST  was  written  primarily 
to  generate  histograms  and  histogram  statistics  for  the  clutter  data. 
The  basic  tasks  of  the  program  were  as  follows: 

a)  Read  and  decode  a data  reduction  tape 

b)  Form  summation  of  time  histograms  and  normalized 
histograms  for  forming  range  gate  and  total  histograms 

c)  Extract  pertinent  statistical  parameters  for  the  time 
histogram  (timehist) , range  gate  histograms  and  the 
total  histograms 

d)  Print  out  statistical  parameters  and  certain  histogreuns 


UNCLASSIFIED 
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F 


e)  Write  data  tape  in  CDC  format  for  further  data  analysis 

® of  histogram  results 

f)  Extract  average  power  per  frame  from  data  tape  and 
write  on  tape  in  CDC  format 

The  program  has  the  following  general  annotated  block  diagram: 


Progrcun  begins  by  clearing  histogram 
arrays  and  read-in  pertinent  arrays. 


It  now  sets  up  for  random  access 
storage  of  the  time  histograms 
(timhist)  arrays  for  each  frame, 
zeros  are  written  into  timhist 
arrays , 


The  program  now  reads  from  the  data 
tape  the  title  and  writes  this 
information  on  tapes  receiving 
histogram  data  and  statistics 


20  is  the  return  point  for  the  next 
timhist  after  the  previous  timhist 
per  frame  and  the  necessary  data 
has  been  extracted  and  recorded. 


Reads  in  the  first  block  and  extracts 
range  gate,  histogram  frame  number 
and  altitude  from  the  first  word 
in  record.  An  end  of  data  causes 
the  program  to  jump  to  80. 


Reads  from  mass  storage  a new 
timhist  array. 
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When  new  block  contains  a change 
in  range  gate  (r.g.)  tag,  the 
program  jumps  to  60  to  output 
range  gate  sum  histogreuns  (R.G.  hists) 
and  statistics. 


The  progreun  normalizes  per  logic 
control  given  at  input  and  either 
normalizes  the  time  histograms, 
R.G.  hists,  and  total  histograms 
(Tothist)  or  only  the  Tothist. 


The  progreun  now  reads  from  the 
data  tape  and  extracts  one  word 
giving  the  average  power  for  the 
previous  frame. 


The  average  power  is  printed  out 
onto  the  lineprinter. 


The  progreun  now  extracts  600  words. 
These  are  the  averages  over  600  FFT 
frames  contained  in  the  previous 
histogram. 


The  600  FFT  averages  are  now  put  on 
MIK  TAPE  and  also  put  out  on  line 
printer. 


The  timhist  array  of  the  previous 
frame  is  written  to  mass  storage. 


32  wqrds  are  now  extracted.  These 
are  the  average  power  for  the  32 
doppler  cells  over  600  FFT' frames 
in  the  previous  histogreun  frame 
record. 
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The  32  doppler  cell  averages  are 
printed  on  the  line  printer. 


The  32  doppler  cell  averages  are 
written  on  MIK  TAPE. 


There  are  logic  switches  which 
allow  selected  histograms  to  be 
out  on  line  printer  and  out  on 
MIK  TAPE 


When  the  end  of  information  is 
encountered  in  a new  block  of  data 
the  program  branches  to  this  point 
to  put  out  last  range  sum 
histograms  data  and  statistics. 

The  time  frame  sum  histograms  and 
the  total  sum  histogram  with 
statistics  are  also  printed. 


The  program  now  reads  from  mass 
storage  the  time  frame  histogram 
data. 


The  array  of  frame  averages  are 
written. 


The  information  requested  by 
control  cards  are  written  out 
on  tape. 


/ 


1.1.2  Probability  Plots  (PROBP) 

Probability  plots  are  generated  from  a data  tape 
that  is  written  by  the  histogram  program  (TAGHIST) . 
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The  probability  density,  distribution  and  mean 
value  for  each  range  gate  are  read  sequentially  from  the  data  tape. 
After  the  data  for  each  range  gate  is  read,  two  running  sums  are 
accumulated.  The  first  is  a sum  of  normalized  range  gate  histo- 
grams, normalized  range  gate  by  range  gate,  and  the  second  is 
just  a sum  of  the  range  gate  histograms.  When  the  data  for  the 
last  range  gate  is  processed,  the  probability  density  and 
distribution  of  the  total  histogram  is  read.  This  data  is  generated 
in  the  histogram  program  and  is  already  normalized  histogram 
frame  by  histogram  frame. 

"TOTAL-N"  type  plots  are  generated  from  this 
histogram  normalized  data.  The  accumulated  range  gate  "normalized" 
data  is  then  processed  to  produce  "TOTAL-A"  type  plots  and  the 
accumulated  "Un"-normalized  range  gate  data  is  processed  to  pro- 
duce "TOTAL"  type  plots. 

1.1.3  Statistics  Plots  and  Processing  (STATP) 

Statistics  plots  and  the  associated  processing 
is  done  from  a data  tape  written  by  the  TAGHIST  program. 

The  necessary  statistics  are  read  from  the  data 
tape  for  each  frame  and  each  range  gate  and  are  stored.  The 
stored  data  is  scanned  to  eliminate  frames  and  range  gates  with 
incomplete  sets  of  statistics. 

Plots  are  generated  of  the  statistics-vs-f rame 
number  and  statistics-vs-range  gate.  The  mean  is  removed  from 
each  data  set  and  a normalized  set  of  plots  is  also  generated. 

Statistics  of  the  mean  are  computed  for  the 
data  set  where  the  mean  was  computed  in  each  histogram  frame. 
Statistics  are  also  computed  for  the  mean  normalized  median, 

10  10  10  ^ and  10  ^ probability  levels.  These  "statistics" 

of  the  statistics  are  then  output  to  the  printer. 
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1.1.4  Average  Plots  (MEANR) 

Average  plots  are  generated  from  a data  tape 
^that  is  written  by  the  histogram  program.  The  histogram  program 
reformats  this  data  so  it  is  compatible  with  the  CDC  computer 
system. 

Data  is  read  from  the  data  tape  a range  gate 
at  a time.  One  hundred  FFT  frames  are  averaged  at  a time  to 
produce  an  output  data  point.  At  the  completion  of  a range  gate, 
the  next  set  of  range  gate  data  is  read. 

The  average  range  gate  data  is  plotted  a range 
gate  at  a time.  Each  range  gate  is  offset  from  the  preceeding 
range  gate  by  2dB  so  the  average  plots  will  not  overlap. 

1. 2 Hit  Software  Flow  Charts 

Overall  software  flow  for  hit  processing  is  summarized 
in  Figure  D-2.  The  DT-1  hit  tape  received  from  the  Data  Reduction 
Facility  is  sorted  from  sequential  range  gate  to  sequential  time 
by  use  of  mass  storage  devices  (i.e. , computer  disk  storage) . The 
sorted  hit  tape  is  used  for  all  subsequent  software  flow. 

The  three  main  paths  indicated  summarize  software  flow 
for  Hit  Map,  Hit  Counts  vs.  Time  and  Conditional  Probability 
Map  Processing.  A data  processing  program  and  plot  program  are 
shown  within  each  path.  Various  summary  printouts  are  generated 
for  use  as  diagnostic  tools  at  each  step  in  processing.  In 
addition,  optional  printouts  and  special  dump  options  are  built 
into  the  software  for  special  analysis. 

Detailed  flow  charts  for  Hit  Map  processing  in  Figures 
D-3  and  D-4  describe  input/output  operations,  processing  algorithms 
and  plot  steps.  Typical  hit  map  printout  is  presented  in  Figure  D-5. 
Hits  in  each  100  ft.  X 100  ft.  cell  are  listed  along  with  sums 
over  range,  time  and  a total. 
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Figure  D-3  Hit  Data  Processing  for  Hit 
Map  and  Large  Hit  Detection 
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Figure  D-5  Hit  Map  Line  Printer  Output 
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Hit  Counts  vs.  Tine  data  processing  is  indicated  in 
Figure  D-6.  Plot  program  flow  is  detailed  in  Figure  D-7 . Hit 
Counts  in  each  FFT  tine  fraune  and  range  gate  are  written  to  an 
output  tape.  The  plot  progreun  contains  processing  software 
sufficient  for  plotting  fine  and  coarse  Hit  Counts  vs.  Time  from 
the  same  data  tape. 

Processing  for  generation  of  Conditional  Probability 
Maps  is  detailed  in  Figure  D-8.  Data  processing  for  construction 
of  the  Condiyional  Probability  Cube  is  defined  in  the  left  half 
of  the  flow  chart.  Plot  processing  is  indicated  on  the  right 
half  of  the  flow  charts.  Options  are  built  in  for  plotting 
cuts  by  time,  range  and  doppler  on  all  cuts.  The  time  collapsed 
array  and  normalization  array  are  standard  output  for  any  choice 
of  cuts  through  the  Conditional  Probability  Cube. 
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Figure  D-6  Hit  Counts  Versus  Time  Data  Processing 
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Figure  D-7  Hit  CoVmts  vs  Time  Plot  Program 
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Figure  D-8  Software  Processing  for  Conditional  Probability  Maps 
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APPENDIX  E 
SIMULATION  SOFTWARE 

1 . 1 simulation  Structure 

The  simulation  of  the  TAGSEA  system,  clutter  model,  and 
analog  data  reduction  was  designed  to  provide  accurate  results 
while  minimizing  the  eimount  of  complex  calculations  required. 
Considerable  effort  was  devoted  to  improving  efficiency  since 
literally  millions  of  numbers  need  to  be  generated  to  simulate 
a run.  The  simulation  was  specifically  designed  to  model  the 
TAGSEA  system  only.  Ideas  and  sections  of  previously  existing 
clutter  simulations  were  utilized,  however,  and  the  simulation 
was  configured  in  such  a manner  that  generalization  and/or 
modification  to  another  system  could  be  accomplished  quite 
easily. 

A flow  chart  of  the  simulation  is  shown  in  Figure  E-1. 
The  program  starts  by  reading  a set  of  inputs  which  define  the 
geometry  and  various  system  parcuneters.  Then  various  system 
parameters  which  can  be  considered  "stationary"  over  a short 
period  of  time  are  generated  and  placed  in  arrays.  The  clutter 
model  is  then  used  to  provide  values  with  the  correct  amplitude 
and  distribution.  These  values  are  then  passed  through  a model 
of  data  reduction  analog  processing  which  modifies  the  data  on 
the  basis  of  the  best  estimate  of  the  clutter-independent 
; pareuneters.  These  clutter  values  are  then  written  to  two  tapes 

which  are  equivalent  to  the  digital  tapes  normally  provided 
by  the  data  reduction  facility.  After  the  values  have  been 
written  to  the  tape  a decision  is  made  to  stop  of  not.  If  the 
run  continues  a further  decision  is  made  as  to  whether  to  update 
i the  "stationary"  pareuneters  or  not.  At  the  end  of  the  run  the 

I tapes  are  sorted  to  the  correct  order  and  passed  on  for  process- 

ing by  the  standard  TAGSEA  software. 
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Figure  E-1  TAGSEA  Simulation  Flow  Chart 
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The  transfer  function  block  diagram  for  the  simulator 
is  shown  in  Figure  E-2.  Shown  in  this  diagreim  are  the  effects 
modeled  in  the  simulation  and  the  neeunes  of  the  subroutines 
(discussed  in  the  next  section)  where  the  modeling  takes  place. 

1. 2 Detailed  Description 

The  following  subsections  contain  the  detailed  code  and 
description  of  the  TAGSEA  simulation.  A listing  of  the  main 
program  is  provided.  This  listing  contains  a definition  of  the 
program  inputs  and  outputs  and  a key  variable  dictionary.  It 
also  defines  the  overall  structure.  A description  of  each  sub- 
routine is  also  included,  giving  the  inputs,  outputs,  and 
modeling  assumptions  used  in  each  routine. 

1.2.1  Progrctm  Listing 

Table  E-1  (in  four  parts)  is  a reproduction  of 
the  simulation  program  printout  including  comment  cards.  It  is 
the  basic  function  of  this  program  to  implement  the  flow  chart 
previously  described  and  to  call  up  the  subroutines  which  are 
the  subject  of  the  succeeding  subsections. 
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Figure  E-2  TAGSEA  Simulation  Transfer  Function  Diagram 
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MlTtHlSTO.TAPCl ItMfcOS  tTAPE98 
iTAREl0>HlT.TAPEI2>HISToiTAPE9«;M60S> 
XPOS(3?*l6)t  YP0S(4Ptl6) 


H1sTARV(166o)« 

FHmPhR(600,16>. 


PLvCi./4A0«/', 

TPLALTO/10oo*0/i 
■DOIr/90.0/, 
ONaNoe/IOO.q/* 
0OpOEF/4000.0/« 
OTmIsTO/5,0/  I 


plaltOi  plvelp* 
MUOlpt  LiR* 

RAnOE#  OOPOPP» 

s»tfNeot 


»TOR«UOfUi 

p0PAV0(32,i4) 


PLALTO/0.0/  ■ 
rPLVbLO/lOoe.O/ 

LAM/U.I/  

UOOP/i09*0/ 

foF/i.O/' 

SHFptQ/141no,0/ 


TPLALTO 

pRF 

ToF 


phospam  to  Simulate  taosea  oaja  oeneration 


PRC OR Am  inputs 

PLALT  altitude  oF  plane  «fT» 

PLVEL  YtuOclTT  oF  PLANE  (FT/SECJ  

PLALTO  rate  CF  ChANgE  OF  PlANE  AltITuOE  CFT/SEU) 

plvElo  rate  of  change  of  Plane  vflOCitv  (Ff/SEU2» 

tplaltd  period  of  Plane  altitude  variation' (SeC) ' 

TPLYELo  PtRloO  Of  PL*ne‘ velocity  vaRXATION  (SCC» 

NOVEL  aInO  velocity  (FT/^eC) 

■ooiR  aino  direction  (cluckmIse  fROm  Plane  velocity')' (ocal 

LAM  AAvElEnOtm  (FT) 

PRF  transmitter  phf  (mZ  ) 

ORANGE  RAn(»E  rate  MIOTM  (FT)  

OOOP  FFt  cell  mIOTH  (HjJ 

RAKOE  MINIMUM  gated  RANqE  IFT> 

OOPOFF  IF  OoPPLcR  OFFSET  ImZ»  - 

TOF  time  of  FlIOmT  (SEC) 

otgeo  time  Increment  to  OpoaTe  gfomctry  isec» 

otmisto  time  Increment  to  nqite  histogram  RlcoruS  (Sec) 

ShFREU  SmmPlE  And  hold  FREOUCNCY  (MZ) 

MITTHR  MlT  THRESHOLD 


UNCLAsIlFIED 


O O O O O O U U o o o o w u o o i>uo  o o o ouo  uuouu 


UNCLASSIFIED 


TABLE  E-1 

TAGSEA  SIMULATION  PROGRAM  (cont'd) 


PXOgRAH  r«0StM  T4/T4  tipT.I  FTN  4.b«420  76/( 


PROliRAM  o^TPUTb 

file  MfSTO  H1sTo6R*m  HEcORDi  FriR  TAOSf*  PROCESSINO  — ' 

file  MtT  sorted  hit  HECOROS  fop  TAfisEA  PROCESSINU 

major  internal  variables 

CELLA  U,  16  array  USEU  Tq  STORE  SEr I-STA T IONARy  EFPECTS 

CELLAB  JR,  16  array  CONTaININO  TOTAL  AND  PHOCESSED  CLUTTCR  PETI 

vEL  total  pLanE  velocity  — 

onnor  nohmalNinO  nAiN  e<>oivalent  to  calibration  Si^iial 


InITILIZaTION 

constants  and  initial  values 


PIs4,O#ATAM1«0|  _ _ 

T-0.0  " ■■  ■ ■ ■ 

TB0HI»2.64PI 

TL»ST»o.o 

initialize  files 

REpIND  is  ■■■  - ■ ' “ 

CALL  OpEKMS(lI<INnEX<B0l'0> 

REPIF^  12_  _ _ 

REaO  inputs 

READIS.IM  ~ 

WRITE<6iIN) 

define  BeOPETHY  

continue 

ALT"PLALT»PL*LTd*sIN»TvOPt4T/TpLalT01 
VELP«PLVEL«PLVtLO»SlN(TPOPI«T/T^LV£LOI 

VELSO»vELP«VtLP**OVELApOVEL-2.0*VELP*'<nvEL«COS«  URO»*HOOIR'| /TWOPl) 
VEL"SOQT(VELSa)  * ■" 

clear  aRravs  

00  10  I*lil* 

DO  10  jai*32 

CELL»(J*I)«1.0  ' ~ 

|0  CONTINUE 

ONNORPItO  _ . 

FIND  area  cf  each  Patch 
AND  X anc  t positiOh  for  lrte*  Use 

AND  mOdIFT  cell  array  By  

•area 


UNCLASSIFIED 


UNCLASSIFIED 


TABLE  E-1 

TAGSEA  SIMULATION  PROGRAM  (cont'd) 


tAOsIM  OpT«l  FTN  4.6*420  76/0 


C 

CALL  TrU«RE  A ( CtLL« t XPOS*  VPOSt  A|.T, R«N6E.  vLL. LAM.  DUOP. URANGE, rNNOR) 
C 

C FIM)  R«»*  FOR  LACh  rIng  “ 

c ARC  MODIFY  cell  Array  by 

C /M««4  _ _ 

c 

c 

call  TrUR4(CElLA.RAn0e,0R6N0E.0RN0R) 

C FIM)  ANTfNAAGAIN  FO»  E*Cm  RATcH  ' ~ ' --  

C modify  cell  ARRaY  Bv 

C mGaIn 

C 

CALL  TRUANYlCELLA.XPOs.YPoS.ALTtPLVEL.WoVEL.MOOIri.VCL.GNNOR', 

C 

c calculate  SI8MA  zero  AS  * functIom  

c cF  graziko  angle 

C MOOlrY  CELL  array  By  _ _ _ 

c MSitto  ‘ ' ' 

c 

call  TRUSlGO»CtLLA.ALT,RANOEil)RANOE»6NNoR) 

C 

C calculate  sin  A/  X FOR  TMe  AMPi.to  SPECTuUM 

C modify  cell  array  By  _ . . ^ 

C MSlNx/X 

C 

CALL  TRUSlAX(CtuLA.PRF,OOoP,OOpOFF»  _ _ _ _ 

C . — . 

C AT  This  point  THb  CELL  ARRaY  hAS  all  TMt 

c semi-statIcnary  Parameters  accoOrteo  For  _ ' 

C Nt*T  IS  THE  generation  of  The  RArOOM  PARAMETERS 

C AND  The  BEMOYEL  Of  TElE  EsTiMATfcS  OF  TerM* 

C 

C - - 

300  CONTINUE 

C _ _ 

c generate  Random  amplItuOCs 

C PLACE  INTO  CELlAB  ARRAv 

C CElLaBbCElLArRANDoM 

c 

CALL  GeNaMP (CELLA.CELlaB.XPOSi Yfos. T.VEl> 
c 

c simulate  Tele  bat  processing 

c 

c since  this  is  Done  repEataTalT 

C IT  MILL  BE  oOnE  In  OnE  PASS 

c 

c „ _ 

Call  RmVtEL (CELLAR. aLt.OOPOFF.gNnOR) 

C 

C DATA  Is  NON  In  a form  SUITABLE  FpR  OUTpgT  TO  TAPES 

c simular  to  tele  TaPes  From  bat  ‘ 

c check  the  uaTa  and  nrITe  the  Tapes 


UNCLASSIFIED 


UNCLASSIFIED 


TABLE  E-1 

TAGSEA  SIMULATION  PROGRAM  (cont'd) 


PH06RAM  TAOSIH  T4/T«  OPTsT 


FTN  «,b*«20  Tb/0 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


FlfibT  (jUdNTIie  TO  1o24  AMPLITUDES 

CAUL  NoRMlCELLABtNCEULABiMAXtMiN, — - ~ 

"RITE  SORTED  rtlT  TAPjf __  

CALL  "rTpIT (CELlAR«n1TTMR,MAX»m1N) 

calculate  and  "Rite  MjsTO  tape  ' 

CALL  "RThIST  (CELLARiOTHlS!OtMA*,HIN,ALTlPRAME)  

PHOCESSInO  pom  THIS  ObeUL  is  no"  FINISmcO 

CHECK  rOR  END  OF  RUN  anO  L^OP  a*CK  TO  START  ANOTHER  UMELL 


if»t,ot,tofi8o  to  Is 
T"T*U8./ShFHEO 

CHECK  ir  time  To  uPd*te  oeohetrt 

IF(T»TlAST,lT.0t6C0»  00  To  300' ■ 
tlast-t  . 

00  T6  

END  OF  FlIOHT 

continue  : ■ ■ 

CLOSE  mIsTC  «|T/IlES_ 
ENOPILe  10 

ENCPILe  10 

HISTO  tape  must  "e  sOrtEO  To 

ranoe’TIme  prom 
TImC.RaNoe 


call  SoRtHIS(FHaHe) 

ENOPILe  12 
ENOPlLe  1* 

END  OF  Run 

STOP  1 

End 


UNCLASSIFIED 


UNCLASSIFIED 


1.2.2  True  Area  Subroutine 

The  true  area  subroutine  of  Table  E-2  responds 
to  the  call  up  TRUAREA  of  Table  E-1.  The  area  is  computed 
using  a trapezoidal  approximation  of  the  surface  range-doppler 
cell.  All  four  corner  locations  are  computed  and  the  center 
found  by  averaging  the  corners. 


TABLE  E-2 

TRUE  AREA  SUBROUTINE 


sutiinuTiNe  r*u«ac«  7«/(«  o»t«i  ftn  *.6«4?o  t«/o 


subroutine  TRUiRCXcCLLRtyROSiVBOSttLTiMRNSEfVFUtLaMiOOORtORANOet 
f BNNOR) 

c 

C THIS  SURROUTINC  rjNOS  Tit  TRIlf  *RE*  0*  tME  aROUNO 

C CElLS  «LSa  FINOS  The  I ANS  V POSITION  0»  the  center  of  the  cells 

c 

C NOaIFIeS  sell*  or 

c 

c 

C NOoIFlrO  RARRNETEWS 

3?. I*  a*ouno  cell  array 

3?.  16  » POMtION  of  center  Or  CE|  L 
3P.16  vROSITtON  OF  center  «F  CELl 
N0R“Aui71Nb  aAIN  for  tele  ►'RoCESsINO 

inputs  PEOUfiED  - 
ALTITinP  or  dLANEIFTI 
nInINum  RANSr  RECEIVED (FTI 
total  vELOtlTV(FT/SCC» 

VAvELEN''H»rT» 

OFL^A  OOPPLE  between  cells «M*> 
delta  pANOt  between  cells (Er, 

c 


C CEi  LA 

C XPoS 

C YRoS 

C ONnOR 

c 

c 

C ALT 

C RAnOE 

C VEi' 

C LAm 

C OOoR 

C ORtNOE 


E-9 


UNCLASSIFIED 


UNCLASSIFIED 


4 

1.2.3  R Subroutine 

As  with  the  previous  subsection  and  those  to 

4 

follow,  the  R subroutine  (shown  in  Table  E-3)  is  called  up 

by  the  program  of  Table  E-1.  It  is  the  function  of  this  particular 

4 

program  to  compute  the  R term  to  the  center  of  each  cell. 


TABLE  E-3 
r'*  SUBROUTINE 


SOBBOUTINC  TBUBi  ?*/I4  0BT«1  - - PYN  «.«*4r0  T*/t 


SUhROUTINB  TByB4JCeLL4i***ii9ttOB»N«t.»»m««S  • 
c 
c 

c THIS  suRboot^cn  riNOs  b*44  Bob  c»cm  eru 

C TO  THt  CCBTCB  OB  T«t  CttC 

C 

C MOnirifO  B»B»wfREBS  - . 

C cell*  WORKINQ  ORUUmD  CELL  ARRAY 

C ONmOR  NORMAlTTINO  nAlN  BqR  REMOVlNo  TC|.E  CONSTANTS 

c 

C INPUTS  rcbuireo 

C RAmOE  RANOE  TO  EUOp  OF  FIRST  CCLW(fT| 

C 0R»N8E  delta  RANOt  pEThcEn  CELLS  itT* 

C 

c 


E-10 

UNCLASSIFIED 


UNCLASSIFIED 


1.2.2  True  Area  Subroutine 


The  true  area  subroutine  of  Table  E-2  responds 
to  the  call  up  TRUAREA  of  Table  E-1.  The  area  is  computed 
using  a trapezoidal  approximation  of  the  surface  range-doppler 
cell.  All  four  corner  locations  are  computed  and  the  center 
found  by  averaging  the  corners. 


TABLE  E-2 

TRUE  AREA  SUBROUTINE 


SUtOnuriNE  TRU«aC«  T4/t*  0»T«|  rTN  4.6*4;0  T6/0 

SUnROUTlN*  rRUt4C4tcCLL4fv40^tV»OS«4LTfM*MAEtvFLtL«Mi0004,DR*NOe* 
f ANNORI 

c 

C tMIS  SURROUTIWE  RT*'t)S  TiE  TRHf  4RC4  0'  tME  aROUND 

C CElLS  «LSa  riNAS  The  X 4Na  Y position  0^  the  center  or  the  cells 

c 

C MOniriES  SELL*  BY  a*RE» 

e 

c 

C MOAirtrO  R4R4METERS 

C CEiL* 

C XPoS 

C YPaS 

C ONigOR 

C 

e 

C »Lt 

C RAnOE 

C VEi’ 

C L*m 

C OOaP 

C ORiNSE 

C 

c 


! 


E-9 

UNCLASSIFIED 


3?. 14  aPOUNO  cell  4RP4T 
3?. 16  X position  or  CENTER  Or  CE^L 
3P.14  vPOSITtON  nr  CENTER  wr  CEL^ 
N0P“4ui?INto  a*IN  roR  TELE  rRnCESsINO 

INPUTS  REOUfiEO 
ALTiTunP  or  ALANE(rT) 
nINI**um  R4N6r  RECEIVEO<rT» 

Total  vELOtlTYirT/SECi 
P4YCLEN'*HfrT» 

OEl’4  oopplE  between  cells (Hr) 
delta  RANOt  aETWEEN  cells  irr, 


UNCLASSIFIED 


1.2.4  Antenna  Subroutine 


Antenna  gain  is  found  to  the  center  of  each 
patch  using  an  analytic  formula  for  the  antenna  gain 

G - U - 


11  = peak  gain  in  dB 

~ 4.61  X 10~^  5eg2  azimuth  rolloff  constant 

K„.,  = 4.70  X 10  -3 — 2 elevation  rolloff  constant 

EL  deg^ 

Gain  is  assumed  constant  over  the  patch.  (See  Table  E-4) . 


TABLE  E-4 

ANTENNA  SUBROUTINE 


SUMOUTlNt  TI*0*nT  ■ 0^T»| 


fTN 


SUpOOurlNB  TRUANT  (CELL*  tXPOS  ,VP0St«LTi''LvCL8WnvCLiMD0IRtVEL»0NN0R^ 
P I 

THIS  SuBROWTtNE  CALCULATES  ThE  AnTENN*  RAIN  ' 

TO  THE  CENTER  ftP  E»cH  RATfH 

TAPES  into  account  THE  CRJS  anOLE  CAUSEU  hV  ThE  MInO 


MOpIEIEO 

CElLA 

SNdiOR 


PARKINS  ARNAy  • SAIN 

normalizing  sAlN  USED  EOR  TEL'E  Bir  REMOVAL 


INPUTS  RESUIREO  

XPoS  3?tl6  aRRAT  nP  * CENTERS  Ot  pATCmESIPTI 

VPnS  3P«16  aRRAT  ftW  V CENTERS  Ot  PATCHES IPT) 

ALT  ALTITueE  op  pLANEIPT) 

PLVEL  velocity  op  tHE  AlRPLANEfET/SECr ' 

VOvEL  PINO  velocity CET/SECI 

HOnIR  PINO  OiRECTInNIOEOi  CLOCKHlSp  ERoM  PLANE  VELOCITV 

VECTaR  - - - 

VEl  total  velocity iet/seci 


internal  tONSTANTS 

KAT  ATMUTh  ANTINnA  ROOLOEE  COEt 
KEl  elevation  AN^  ROOL  OEE  COEP 
OAnS  ant  OERRCSPI(bN  ANSlC 


E-11 

UNCLASSIFIED 


UNCLASSIFIED 


1 


1.2.5  Mean  Backscatter  Subroutine 

Mean  backscatter,  , is  calculated  as  a 
function  of  altitude,  range  and  doppler  cell  position  for  each 
surface  cell  using  the  program  constants  which  were  inputted 
to  the  overall  program  (see  Table  E-5) . 


TABLE  E-5 

MEAN  BACKSCATTER  SUBROUTINE 


sowtoortMC  T»u8f«o  - t>rr«i --  ■ rtN  4,6**80  - 


SUllWOUTlNl  T»OSieo  fceLL*tSLTt<»»N«etOS*N*€;eH«8«t 

c 

C THtS  SuBRBUTInC  C4LC'JL4Tf,  STQM*  ?ERO  i FUnCTION  OF 
C »«i«OE  OiTE  --  - 

C 

C ROniriEO  BARAMETEftS 

C CClLA  3E,16  wORKlNft  PARAMETER  Cett*  

c 

c 

C INtCRN»L  MNSTANTS  

c SI AO  16  array  0^  aIOMA  EERO  CON»TaNTS 

c 


1.2.6 


UNCLASSIFIED 


1 


Sin  x/x  Subroutine 

Calculation  of  the  transfer  function  of  the 
sample-and-hold  implementation  over  the  16  range  gate  arrays 
is  performed  by  the  subroutine  listed  in  Table  E-6. 


TABLE  E-6 

SIN  x/x  SUBROUTINE 


lUMOOTiNE  TltUSf»<*  FTN  4.4*4S0' 


- sobbootin*  tbusinx  ecELL^f ii»r t000Bt00i*0f?r - 

c 

C TMtS  SUBBBUTiNf  ADDS  IN  TmE  errfCTS  or  »/m  0^fl»*TION 
t ON  THE  OENtdjTEO  SrECTBOT  - 

C MOnlFlEO  IAR4METERS 

C CEtt*  3?»16  wOBKINb  cell  4m»4r  •»tNX/X  

c 

C INbuTS  RCtiUXBEB 

C BRr  SYSTEN  ruriMri  - 

C OOoP  DOPPLER  BETWpBN  FFT  CELLS 

C _OOpOFF  2ER0  0ol*PLtR  OFFSET  IN  FFT«Hj» 


E-13 

UNCLASSIFIED 


t 


tS/< 


4. 


A 


r 


UNCLASSIFIED 

1.2.7  Gain  Normalizing  Subroutine 

In  the  data  reduction  processing  of  the  collected 
clutter  data,  the  gains  for  the  16  range  gates  and  32  doppler 
cells  are  normalized  to  restrict  the  dynamic  range  required  in 
further  processing.  The  gain  normalizing  subroutine  of  Table  E-7 
performs  the  same  function  for  the  simulation. 


TABLE  E-7 

GAIN  NORMALIZING  SUBROUTINE 


UNCLASSIFIED 


1.2.8  Integer  Array  Subroutine 

Clutter  power  over  the  array  of  cells  in 
essentially  continuous  (analog)  form  is  normalized  and  arrayed 
into  1024  bins  corresponding  to  the  like  bins  of  the  actual 
clutter  data  histogreuns  (see  Table  E-8)  . 


TABLE  E-8 

INTEGER  ARRAY  SUBROUTINE 


SUBROUTlNe  NORM  ORTal 


rXN  4,6*4»0  Tt/( 


SUBROUTlNft  NO»M(CCLL***'*eRLL«ll*M«KtM*N»  

OlMCNStON  CeUAB  02,161  *NfeLL6R'|3>*l*l 

TMtS  subroutine  CHAnBCS  6 RC»U  ROMEO  ROMRf  IWtO 

*N  INTEOfO  *RR«V 

6.L0MS  OOOVISXON  FOR  6 «»tN  CHANOC  «T  T«E  SR%  TIME 


FINOS  THE  MIN  »N0  TmE  M*«  IN  THE  »00»V  — 

TRUNCATES  THE  AMPLITUDES  tO  MITmIN  ^ TO  I#81  ; 

MOnlFlEO  MARAMETEOS 

NCfLLAR  32,16  iNTEMEo  ROMEO  AOOAV 

lEOUAVALENLEn  TO  CELLAR  IN  OitO  pOOBOiUlT 
MA«  maximum  AMMLtTUOE  in  the  ANRAT 

MIn  minimum  amplitude  in  the  AHOaT 


internal  tONSTANTS 

OAlN  bain  MoOIFlCiTTON  ALLOMED 


i: 

I 


UNCLASSIFIED 


UNCLASSIFIED 


1.2.9  Hit  Subroutine 

Using  the  inputted  threshold,  the  hit  sub- 
routine of  Table  E-9  sorts  the  simulated  clutter  returns  and 
write  all  hits  greater  than  the  threshold  onto  a hit  tape. 


TABLE  E-9 
HIT  SUBROUTINE 


■ugwuTrNt  whtmtT 


nu  4.**4>o 


SUnROUTlNf  WRTHITICtLL*««MlTTH1»,**«*»1ITNt 
dimension  CELLABOZ,!*! 

INtEOER  BELLARtHITTHR 


THtS  SUBRRUTINE  writes  UAqOE  amplitude  NCfURNs  TO 
» sorted  RIT  T*PE  for  further  RROCESSIN® 

BY  normal  TAOSFA  SOFTHANE 

HOoIFIeO  RARARETERS  “ ' “ 

none 


INPUTS  REBOlRtO  “ ' 

CEiLAB  3?#16  iNTCWEo  ARRAY  OF  NOR"AlIZEo  AMPLITUDES 

ha;  LyROEST  AMMLrTUOE  IN  THE  AHR*T,BpFORE  TRUNCTION 

MIn  smallest  AMP^ITUEE  in  the  ARpAYIpEPORE  TRUNCTAION) 

internal  constants 

IFpAME  internal  eft  FRAN  COUNT ' 


E-16 


UNCLASSIFIED 


UNCLASSIFIED 


1.2.10  Histogreon  Subroutine 


Histogreun  tapes  corresponding  to  those  generated 
from  actual  data  and  using  the  same  format  (except  for  range- 
time inversion) , are  generated  by  the  subroutine  of  Table  E-10. 


TABLE  E-10 

HISTOGRAM  SUBROUTINE 


SUbJIOUriNE  VRTHfST  T4/T4  O^Tal 


rrN  4.6*420 


SUnROOTlN*  W»TMlSTlcfLL*B^OTHI8TO,M»X»M*M;*l.T;rR»*tr)-  

THtS  SUBROUTlMr  GEMEOATtS  ARRAYS  AND  RRAAmETErS  REOUIRCO 
FOR  histogram  record  TAbES  ~ 

outputs  tapes  in  The  ORUEr  riMEt  RANOE#  OoRRLE* 

REoUIRrO  tUPOTS  

CElLAB  3?. 16  integer  ARRAY  CONTAINOmG  AMPLttEDEI 


-ootrots- 

tareii 


random  ACE6S  PILE  eONTAlNG  HiSTOsRAM  STATlSTICt 


MIrTARy  If  HADE  UR  qR  TME  ROLLONINS  ftOUKs 


word 

1 

t 

4 TO 
10»7 
10P8 
»0?9  To 
1620 
16?9  TO 
16a0 


RARAMEtER  

range  OATE 

FRAME  COUNT 

altitude  

amplitude  NIrTOGRAM 
average  POREr  for  TME TRAN* 

POKER  IN  each  FFT 
average  POREr  PER  DOPPLER  %t 


TO  save  SNACE  large  amplitude  values  ARA  sTORrA  »**-* 

VtPIARLE  IENOHT  sequent  I it  ARRSVtSAEl 


E-17 

UNCLASSIFIED 


r 


m 


UNCLASSIFIED 

1.2.11  Range-Time  Sort  Subroutine 

The  range-time  sort  subroutine  listed  in 
Table  E-11  is  required  because  the  simulation  writes  records 
in  the  order  of  time,  range  gate.  Normal  processing  expects  to 
see  the  records  in  the  order  range,  time. 


UNCLASSIFIED 


1.2.12  Amplitude  Generator  Subroutine 

This  subroutine,  unlike  those  described  in 
the  previous  subsections,  makes  use  of  previous  written  programs 
and  thus  the  listing  will  not  be  included  here.  The  description 
which  follows  is  intended  to  provide  information  for  a more 
complete  understanding  of  the  overall  progreun. 

This  subroutine  is  designed  to  model  the  prob- 
ability density  function  (pdf)  of  the  clutter  process  as  well 
as  certain  statistics  of  the  map  means.  The  process  is  modeled 
on  the  product  of  two  processes:  and  X^.  Xj^  is  a slowly 

varying  process  with  a time  constant  of  a few  seconds  while  X^ 
is  independent  from  one  FFT  look  to  the  next  (and  from  clutter 
cell  to  clutter  cell) . It  is  assumed  that  the  pdf  and  auto- 
correlation of  X,,  and  the  cdf  of  X_  has  been  established  to  model 
specific  clutter  data.  For  example,  in  the  simulation  output 
presented  in  the  report,  a 2.2KFT  altitude  look  at  the  ocean 
in  SS4-5  with  vertical  polarization  was  employed. 

The  simulation  made  extensive  use  of  a pre- 
viously developed  uniform  random  number  generator  (URV)  which 
gives  a range  of  numbers  from  0 to  1.  First,  X^  was  generated 
by  comparing  the  output  of  URV  to  the  amplitude  that  X^  must 
have  so  that  the  cdf  evaluated  at  X^  will  equal  the  number  from 
URV.  Gaussian  random  numbers  were  then  generated  by  using  a 
transformation  on  URV  outputs.  From  the  existing  data  a 
Gaussian  pdf  appears  reasonable  for  X^^  as  does  an  exponential 
fit  to  its  autocorrelation  function.  The  process  can  then  be 
simulated  as  the  output  of  a low  pass  filter  excited  by  "white" 
Gaussian  noise. 
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To  illustrate  the  method  used  to  obtain  Gaussian 
numbers  from  a uniform  number  generator,  we  first  construct  a 
Rayleigh  (voltage)  random  variable  and  then  pick  a random  phase 
from  0 to  2ir: 


u = URV(x)  0 < u < 1 
V = URV(x)  1 < V < 1 
D = (-2  log  u)^/^ 

T = D cos  (2ttV) 


T is  Gaussian  with  zero  mean  and  unit  variance.  The  pareuneter 
X is  simply  used  as  a calling  argument  to  the  systems  URV  number 
generator. 


We  now  examine  the  required  analog  filter  to 
yield  an  exponentail  autocorrelation  function.  If  R(t)  is  the 
autocorrelation  function  desired,  we  note  that 


F 


T i 

c u(T)+e  ° u(-t)  ^ = 


T 

c 


where:  F = Fourier  Transform 

T = the  time  constant  on  the  autocorrelation  function 
c 

S (u)  = Power  density  spectrum 
u(t)  = unit  step  function 

We  also  note  that  white  noise  passing  through  a filter  of  the 

form  H(ja)  = l/(jwt  + 1)  will  result  in  the  seune  spectral  slope. 

c 

The  impulse  response  of  H(ja)  has  the  form: 

I (t)  = e "^c  for  t ^ 0 

This  means  the  autocorrelation  function  and  the  low  pass  filter 
time  constants  are  identical. 


E-20 

UNCLASSIFIED 


UNCLASSIFIED 


The  filter  required  was  modeled  with  the  following 
difference  equation: 

X„(t+1)  = A X^(t)  + B X„(t+1)  + C 
m m N 

where:  X„(t+1)  = New  mean  value  from  filter 

Xjj(t+1)  = New  Gaussian  random  variable  with  mean  of 
1.00  and  standard  deviation  of 

XN 

A = A constant  related  to  the  filter  bandwidth 

B = A constant  picked  to  keep  the  variance  of  X„ 
at  a constant  level  (= 

AM 

C = A constant  to  provide  X„(t+1)  with  a proper  mean 

Parameter  "A"  was  developed  by  using  the  impulse  invariant  digital 

filter  approach,  so  if  the  simulation  sample  time  is  At,  then  the 

value  of  A becomes:  A = e . It  can  be  further  shown  if 

B = (1-A^)^'^^  and  C = 1-A-B,  the  X„  will  be  Gaussian  with  a mean 

of  1.00,  a standard  deviation  of  the  driving  function  and 

will  have  an  exponential  autocorrelation  function  with  time 

constant  T . 

c 

In  the  model,  the  value  of  X„  so  obtained  was 

used  for  every  cell  in  the  FFT,  while  the  value  of  Xq  was  allowed 

to  vary  from  cell  to  cell.  The  resulting  amplitude  from  each 

cell  is  "X"  where  X = X„X_. 

M y 
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This  appendix  presents  the  results  of  the  simulation 
described  in  the  previous  appendix.  Throughout  this  appendix, 
the  simulation  output  data  is  compared  to  one  of  two  data- 
gathering  missions  (run  5 flight  6 and  run  4 flight  6)  both  of 
which  are  high  sea  state  (SS-4  or  SS-5)  with  a 2200  foot  air- 
craft altitude  with  vertical  polarization.  Essentially,  four 
different  types  of  simulation  output  are  provided  and  compared 
to  the  flight  data: 

1)  Various  functions  of  the  PDF  of  the  overall  data 

2)  Map  mean  PDF  and  autocorrelation 

3)  Hit  counts  and  hit  maps 

4)  Conditional  probability  maps 

(For  explanations  of  the  terminology  used  throughout  this  sec- 
tion, the  reader  is  referred  to  Volume  II,  Section  9) . 

Beginning  with  the  first  type  of  simulation  output, 
several  measures  of  the  density  function  are  used  to  provide 
the  detail  necessary  for  interpretation  of  the  PDF.  The  following 
figures  are  presented  to  obtain  the  derived  clarity  as  well  as 
to  provide  for  a comparison  to  actual  data: 

Figure  F-1  Simulation-Log  Density  of  TOTAL 
F-2  Simulation-Log  Q of  TOTAL (A) 

F-3  Simulation-Log  Q of  TOTAL (N) 

F-4  Simulation-Log  Q of  TOTAL 
F-5  Run  5 Flight  6-Log  Q of  TOTAL 
F-6  Simulation-Weibull  of  TOTAL 
F-7  Run  5 Flight  6-Weibull  of  TOTAL 
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Figure  F-1  Simulation-Log  Density  of  TOTAL 
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Figure  F-2  Simulation-Log  Q of  TOTAL (A) 
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Figure  F-4  Simulation-Log  Q of  TOTAL 
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Figure  F-5  Run. 5 Flight  6-Log  Q of  TOTAL 
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Figure  F-8  Simulated  Probability  Density  Function 
of  the  Map  Mean 

F-9  Mean  Analysis  Histogram,  Run  5 Flight  6 

F-10  Simulated  Autocorrelation  Function  of 
the  Map  Mean 

F-11  Autocorrelation  Function  of  the  Map  Mean, 
Run  5 Flight  6 


F-12  Mean  Analysis-Simulated  Power  Spectral  Density 

F-13  Mean  Analysis-Power  Spectral  Density,  Run  5 
Flight  6 


F-14  Simulation  Mean  vs.  FFT  Number 

F-15  Mean  vs.  FFT  Number,  Run  4 Flight  16 

F-16  Simulator  Hit  Counts  vs.  Time  (Coarse) 

F-17  Hit  Counts  vs.  Time  (coarse) , Run  4 Flight  16 
F-18  Simulation  Nit  Count  vs.  Time  (Fine) 

F-19  Hit  Counts  vs.  Time  (Fine) , Run  4 Flight  16 
F-20  Simulation  Hit  Map  vs.  Ground  Position 
F-21  Hit  Map  vs.  Ground  Position,  Run  4 Flight  16 
Comparing  the  Log  Q plots  between  the  simulated  and  real  data, 
it  is  seen  that  less  than  a 0.2dB  difference  exists  between  the 
two  thus  indicating  that  the  random  variables  generated  represent 
a close  correspondence  to  the  observed  data.  It  should  be 
emphasized  that  no  information  of  the  time  varying  characteristics 
of  the  process  involved  is  contained  in  the  PDF  plots.  For  this, 
we  must  turn  to  an  autocorrelation  function. 


In  the  mean  analysis  that  was  done,  the  32  X 16  radar 
map  was  developed  for  each  FFT  loo)c  and  the  mean  value  of  the 
resulting  512  points  was  obtained.  This  map  mean  varies  with 
time  and  one  can  describe  the  variation  first  by  a PDF  of  the 
map  mean  and  secondly,  by  a time  autocorrelation  function  of  the 
observed  mean.  For  simulation  purposes,  a Gaussian  fit  to  the 
PDF  was  made,  while  an  exponential  fit  to  the  autocorrelation 
function  seemed  appropriate.  As  explained  in  Section  10  of 
Volume  II,  the  only  characteristic  of  the  mean  PDF  that  was 
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Figure  F-8  Simulated  Probability  Density  Function  of  the  Map  Mean 
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Figure  F-10  simulated  Autocorrelation 
Function  of  the  Map  Mean 
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Figure  F-11  Autocorrelation  Function  of 
the  Map  Mean,  Run  5 Flight  6 
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Figure  F-12  Mean  Analysis-Simulated  Power  Spectral  Density 
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Figure  F-13  Mean  Analysis  - Power  Spectral 
Density,  Run  5 Flight  6 
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Figure  F-15  Mean  vs  FFT  Number,  Run  4 Flight  16 
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Figure  F-17  Hit  Counts  vs  Time  (Coarse) , Run  4 Flight  16 
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Figure  F-21  Hit  Map  vs  Ground  Position,  Run  4 Flight  16 
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maintained  was  the  ratio  of  the  standard  deviation  to  the  mean. 
This  can  be  observed  from  Figures  F-8  and  F-9.  An  analysis  of 
the  autocorrelation  of  the  mean  indicates  approximately  the  scime 
time  constant  for  the  simulation  data  and  the  actual  data. 

Below  the  region  R(t)  < .4,  the  autocorrelation  function  is 
unreliable  due  to  sample  noise.  Figure  F-10  illustrates  a time 
constant  of  291  FFTs  (29.1  cells  on  the  graph)  which,  at  the 
rate  of  111  FFTs/second,  results  in  a "fresh"  sample  of  the  mean 
every  At  seconds  when  At  is  given  as: 


m-FF^/s  =8-25  se<^°nds 

The  autocorrelation  was  developed  by  examining  10240  FFTs  which 
yields  only  11  good  seunples  of  the  mean.  Figure  F-11  is  the 
comparable  plot  for  run  605.  The  power  spectral  densities  of 
the  means  (i.e. , the  Fourier  Transform  of  the  time  autocorrelation 
functions)  are  also  provided  in  Figures  F-12  and  F-13. 


Moving  on  to  the  third  type  of  output,  we  examine  a 
comparison  of  the  mean  of  100  FFTs  (3200  numbers)  for  the 
simulated  and  real  data  (Figures  F-14  and  F-15) . It  can  be 
observed  that  both  have  roughly  the  same  spectral  constant  and 
amplitude  variation,  but  the  actual  data  decorrelates  somewhat 
farther  from  range  gate  to  range  gate.  Comparing  the  mean  of  the 
range  gate  versus  time  to  the  slope  of  the  coarse  hit  counts 
versus  time  plots  for  the  simulated  data  and  then  for  the  actual 
data  (Figures  F-16  and  F-17) , we  can  note  a correlation  in  both. 

A finer  grain  hit  counts  versus  time  plot  is  also  provided 
(Figures  F-18  and  F-19) , and  again  very  similar  trends  can  be 
observed.  The  hit  map  versus  surface  position  plots  (Figures 
F-20  and  F-21)  both  show  less  correlation  to  the  mean  versus 
time  plots  then  to  the  hit  count  plots. 
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Finally,  the  conditional  probability  map  comparisons 
are  presented  in  the  following  10  graphs: 

Figure  F-22  Simulation-Relative  Time  0 

F-23  Run  5 Flight  6 - Relative  Time  0 
F-24  Simulation-Relative  Range  0 
F-25  Run  5 Flight  6-Relative  Range  0 
F-26  Simulation-Relative  Doppler  0 
F-27  Run  5 Flight  6-Relative  Doppler  0 
F-28  Simulation-Time  Collapsed  Array 
F-29  Run  5 Flight  6-Time  Collapsed  Array 
F-30  Simulation-Normalization  Array 
F-31  Run  5 Flight  6 - Normalization  Array 
Scanning  through  all  these  figures,  we  note  a strong  similarity 
between  the  simulated  and  observed  data.  No  exact  measure  of 
simulation  vs.  actual  runs  is  possible  but  the  texture  of  the 
pairs  of  plots  shows  that  the  simulation  yielded  results  which 
matched  the  form  and  shape  of  the  conditional  probability  maps 
almost  exactly. 

Thus  the  overall  simulation  outputs  have  shown  that  the 
data  gathering,  reduction,  analysis  and  the  simulation  itself 
are  all  verified.  It  shows  further  that  the  clutter  model  used 
is  indeed  a good  model  to  use  for  simulation  which  is  a pertinent 
point  for  use  in  future  missile  system  work. 
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Figure  F-30  Simulation-Normalization  Array 
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Figure  F-31  Run  5 Flight  6-Normalization  Array 


F-40 

UNCLASSIFIED 


